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p I ' Abstract 

We prove global stability of Minkowski space for the Einstein vacuum equations in harmonic 
<**] ■ (wave) coordinate gauge for the set of restricted data coinciding with the Schwarzschild solution in 

the neighborhood of space-like infinity. The result contradicts previous beliefs that wave coordinates 
are " unstable in the large" and provides an alternative approach to the stability problem originally 

1 1 1 solved ( for unrestricted data, in a different gauge and with a precise description of the asymptotic 
behavior at null infinity) by D. Christodoulou and S. Klainerman. 

Using the wave coordinate gauge we recast the Einstein equations as a system of quasilinear wave 
' equations and, in absence of the classical null condition, establish a small data global existence 

— . result. In our previous work we introduced the notion of a weak null condition and showed that 

the Einstein equations in harmonic coordinates satisfy this condition. The result of this paper relies 
on this observation and combines it with the vector field method based on the symmetries of the 
' standard Minkowski space. 

In a forthcoming paper we will address the question of stability of Minkowski space for the Einstein 
vacuum equations in wave coordinates for all " small" asymptotically flat data and the case of the 
Einstein equations coupled to a scalar field. 

S: 

> ! 1 Introduction 

The focus of this paper is the question of global existence and stability for the Einstein vacuum 
equations in "harmonic" (wave coordinate) gauge. The Einstein equations determine a 4-d manifold 
A4 with a Lorentzian metric g with vanishing Ricci curvature 



(N 
> 



We consider the initial value problem: for a given a 3-d manifold £, with a Riemannian metric go, and 
a symmetric two-tensor ko, we want to find a 4-d manifold Ai, with a Lorentzian metric g satisfying 
the Einstein equations, and an imbedding £ C M such that go is the restriction of g to £ and ko is 

*Part of this work was done while H.L. was a Member of the Institute for Advanced Study, Princeton, supported by 
the NSF grant DMS-0111298 to the Institute. H.L. was also partially supported by the NSF Grant DMS-0200226. 

T Part of this work was done while I.R. was a Clay Mathematics Institute Long-Term Prize Fellow. His work was also 
partially supported by the NSF grant DMS-01007791. 
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the second fundamental form of £ in M. The initial value problem is overdetermined which imposes 
compatibility conditions on the initial data: the constraint equations 

Ro - ko) ho* + kol fc j = 0, Vfay - V< fcoj = 0, V» = l,...,3. 

Here -Ro is the scalar curvature of go and V is covariant differentiation with respect to go- The Einstein 
equations are invariant under diffeomorphisms. To have a working formulation one needs to eliminate 
this freedom by fixing a gauge condition or a system of coordinates. 

While the Einstein equations are independent of the choice of a coordinate system, the existence 
of a special or preferred system of coordinates has been a subject of debate |Foj . Historically, the 
first special coordinates were the harmonic coordinates (also referred to as wave coordinates in current 
terminology). These obey the equation D 5 = 0, \x = 0, 1,2,3, where U g = V Q V a is the geometric 
wave operator. Relative to the wave coordinates a Lorentzian metric g satisfies the wave coordinate 
condition if: 1 

(1.1) g ap d p g afl = -s^dpgafl, V/i = 0, ..,3. 

In this system of coordinates, the vacuum Einstein equations take the form of a system of quasilinear 
wave equations 

(1.2) 9 aP d a d pgta/ = Af^(g,dg), Vp,i/ = 0,..,3 

with a nonlinearity Af(u, v) depending quadratically on v. In this particular gauge Choquet-Bruhat 
CB1 was able to establish the existence of a globally hyperbolic development 2 of the Einstein vacuum 
equations starting with an arbitrary set of initial data prescribed on a 3-d space-like hypersurface 
and satisfying the constraint equations. While the result of Choquet-Bruhat and a later result of 
Choquet-Bruhat and Geroch CB-G , establishing the existence of a maximal Cauchy development, 
constructs solutions for any given initial data set, it does not provide any information about the 
geodesic completeness of the obtained solution. In the language of the evolution equations these 
results only show the existence of "local in time" solutions. 

The global results have proved to be by far more resistant. The outstanding global problem, which 
for a long time remained open, and was finally ingeniously solved by Christodoulou and Klainerman 
C-K], was that of the stability of Minkowski space. In simplified language, it is the problem of 
constructing a global solution to the Einstein vacuum equations from the initial data, which is close 
to the Minkowski metric rriavi an d asymptotically approaching the Minkowski space. The initial data 
(£, go,ko) for the problem of stability of Minkowski space is asymptotically flat, i.e., the complement 
of a compact set in £ is diffeomorphic to the complement of a ball in M 3 , and there exists a system 
of coordinates (x%, X2, X3) with r = \J x\ + x\ + x\ such that for all sufficiently large r the metric 3 
goij = (l + 2M/r)5ij +o(r~ 1 ~ <J ), and the second fundamental form ko = o(r~ 2 ~°") for some a > 0. Here 
M is the mass, which by the positive mass theorem is positive unless the data is flat, see Schoen and 
Yau |S-Y| and Witten |Wij . In addition, the data is required to satisfy a global smallness assumption, 
which makes sure that it is sufficiently close to the data (M 3 , 5, 0) for the Minkowski space. 

1 We shall use below the standard convention of summing over repeated indices and the notation d a = d/dx a 
2 For the definitions of global hyperbolicity and maximal Cauchy development see |H-E| . |Wa| 

3 The stability result of |C-K| was proved for strongly asymptotically flat data <?oy = (1 + 2M/r)Sij + o(r _3//2 ), 
fco = o(r- 5 / 2 ). 
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To understand some of the difficulties of the problem we recall that a generic system of quasilinear 
equations 

(1.3) Ufa = M^p^J^^K + cubic terms 

|a|<|/3|<2 

allows solutions with smooth arbitrarily small initial data which blow up in finite time 4 . The key to 
global existence for such equations was the null condition found by Klainerman, |K2j . The small data 
global existence result for the equations satisfying the null condition was established in [HI], 
The null condition manifests itself in special algebraic cancellations in the coefficients Aj^g of the 
quadratic terms of the equation. 5 It can be shown however, that the Einstein vacuum equations in 
wave coordinates do not satisfy the null condition. Moreover, Choquet-Bruhat CB3 showed that even 
without imposing a specific gauge the Einstein equations violate the null condition. 

These considerations led to the suggestion that the wave coordinates are not suitable for proving 
stability of Minkowski space. In fact, considering a second iterate of the equation (|1.2[) . Choquet- 
Bruhat CB2 argued that the Einstein vacuum equations are not stable in wave coordinates near the 
Minkowski solution. All these resulted in the belief that the wave coordinates are unstable in the large 
in the sense that a possible finite time blow up of solutions of the equation (|1.2|) is due to a coordinate 
singularity. 

The global stability of Minkowski space had been proved by Christodoulou and Klainerman |C-K| 
who avoided the use of a preferred system of coordinates and instead relied on the invariant formulation 
of the Einstein equations with the choice of maximal time foliation (or the double null foliation in the 
new proof of Klainerman and Nicolo |K-Nlj ) and utilizing Bianchi identities for the curvature. The 
special structure of the quadratic terms plays a crucial part in the generalized energy estimates which 
form the backbone of the proof but the null condition can not be pointed out precisely. 

A semiglobal stability result was also obtained in the work of Friedrich [FY] . He used the conformal 
method to reduce the global problem to a local one. The approach is invariant and the special structure 
is again exploited implicitly. 

In this paper we revisit the problem of global stability of Minkowski space in wave coordinates. 
More precisely, we consider the data 6 (M 3 , go, ko) with the metric go coinciding with the spatial part of 
the Schwarzschild metric gs = (l + M/r) 4 dx 2 in the region r > 1 >> M, vanishing second fundamental 
form ko for r > 1, and satisfying a global smallness assumption on IR 3 . We prove that for this initial 
data the wave coordinate gauge is stable in the large: the reduced Einstein equations ()1.2|) has a global 
solution g defining a future causally geodesically complete space-time, jH^Ej. The metric g in wave 
coordinates x a , a = 0, .., 3 approaches the Minkowski metric m: sup^g^ \g(t, x) — m\ —* as t — > oo. 

The intuition behind this result is based on the observation that the Einstein vacuum equations in 
wave coordinates p.2|) satisfy the weak null condition. This notion was introduced in |L-R| for general 
quasilinear systems (|1.3|) and requires that the corresponding effective asymptotic system 

(1.4) (d t + d r )(d t -d r )<f> I = r- 1 A i% m {dt ~ d r ) n $j (d t - d r ) m $ K , *j~r0j 

n<m<2 

4 This is in particular true for a semilinear equation \3cf> = (d t 4>) 2 , |J1| . 
5 E.g. D(f> = (d t 4>) 2 — |V;c0| 2 satisfies the null-condition. 

6 The existence of such data is guaranteed by the results of Corvino and Chrusciel-Delay, |Uo| . |C-D| . 
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has global solutions for all small initial data. Here, 

|a|=n,|/3|=m 

The classical null condition states that Aj^ (ui) = and thus implies the weak null condition. The 
asymptotic system (|1.4j) arises as an approximation of l)1.3[) when one neglects the derivatives tangential 
to the outgoing Minkowski light cones, known to have faster decay. The asymptotic equation was 
introduced in [HI] to predict the time of a blow-up for scalar wave equations known to blow up in 
finite time, and was used in jl.2 to find some other scalar wave equations for which the known blow- 
up mechanism was not present. Asymptotic systems played an important role in the analysis of the 
blow-up mechanisms in |Alj . 

In |L-E,j we have shown that the asymptotic system generated by the Einstein equations in wave 
coordinates (|1.2j) has global solutions for all data. In this paper we consider the full nonlinear system 
(|1.2jl . We should note that although the asymptotic system provides useful heuristics about the 
behavior of solutions, in particular the L°° decay of the first derivatives of various components of the 
metric g, it is barely used in our proof of the small data global existence result for the full nonlinear 
equation H1.2|) . While it is tempting to put forward a conjecture that, parallel to the result for the 
classical null condition |C1| . |K2j . the weak null condition guarantees the global existence result for 
small initial data, we can only argue that all known examples seem to confirm it. A simple example of 
an equation satisfying the weak null condition, violating the standard null condition and yet possessing 
global solutions for all data is given by the system 

(1.5) = w • <9 2 ^ + <9V> • D^ = 0, Dw = 

Another example is provided by the equation Ocp = The proof of a small data global existence 

result for this equation is quite involved, |L2j (radial case), [A3] . As we show in this paper the Einstein 
equations (|1.2[) is yet another example. Interestingly enough, at the level of an effective asymptotic 
system the Einstein equations can be modelled by the system (|1.5|) , 

The asymptotic behavior of null components of the Riemann curvature tensor R a ^s of metric g- 
the so called "peeling estimates"- was discussed in the works of Bondi, Sachs and Penrose and becomes 
important in the framework of asymptotically simple space-times (roughly speaking, space-times which 
can be conformally compactified) , see also the paper of Christodoulou |C2| for further discussion of 
such space-times. Global solutions obtained in the work |C-K| were accompanied by very precise 
analysis of its asymptotic behavior although not entirely consistent with peeling estimates. However, 
global solutions obtained by Klainerman-Nicolo |K-Nlj in the exterior 8 stability of Minkowski space 
were shown to possess peeling estimates for special initial data, |K-N2j . 

Our work is less precise about the asymptotic behavior and is focused more on developing a tech- 
nically relatively simple approach allowing us to prove stability of Minkowski space in a physically 
interesting wave coordinate gauge In particular, we rely only on the standard Killing and conformal 
Killing vector fields of Minkowski space and do not construct almost Killing and conformal Killing 
vector fields adapted to the geometry of null cones of the solution g. 

Our proof is based on generalized energy estimates combined with decay estimates. The generalized 
energy estimates are used with Minkowski vector fields {d a ,Q a/ 3 = x a dp — xpd a ,S = x a d a }. For the 

7 For the precise definition see section |S| 

8 Outside of the domain of dependence of a compact set 
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equations satisfying the standard null condition uniform in time bounds on the generalized energies, 
combined with global Sobolev (Klainerman-Sobolev) inequalities, are sufficient to infer small data 
global existence. In our case however the generalized energies slowly grow in time (at the rate of 
t E ) and need to be complemented by independent, not following from the global Sobolev inequalities, 
decay estimates. We derive the latter by direct integration of the equation along the characteristics. 
It is at this point that the intuition from the effective asymptotic system is most useful. We show that 
all components of the metric with exception of one decay at the rate of t . The remaining component 
however decays only as t~ l+£ . Somewhat surprisingly, the glue that holds together such weak decay 
estimates and the generalized energy estimates is the wave coordinate condition 

In this paper we only prove the result for a restricted set of data coinciding with the Schwarzschild 
data outside of the ball of radius one. 9 This allows us to somewhat sidestep the problem of a long 
range effect of a gravitational field. Due to the inward bending of the light rays, solution arising from 
initial data coinciding with the Schwarzschild data outside of the ball of radius one will be equal to 
the Schwarzschild solution in the exterior of the Minkowski cone r = t + 1. 

In our subsequent work we hope to be able to prove the stability of Minkowski space in wave 
coordinates for general data. In addition we hope to show that our method can be also used to treat 
the problem of small data global existence for the Einstein equations coupled to a scalar field. 

Acknowledgments: The authors would like to thank Demetrios Christodoulou and Sergiu Klain- 
erman for their inspiration and encouragement. We particularly benefited from Sergiu Klainerman's 
suggestion to pursue first the problem with restricted data. We would also like to thank Mihalis 
Dafermos and Vince Moncrief for stimulating discussions and useful suggestions. 

2 The main results and the strategy of the proof 

We now formulate the main results of our paper. Our first result is global existence for the Einstein 
vacuum equations in wave coordinates. 

Theorem 2.1. Consider the reduced Einstein vacuum equations 10 

(2.1) n g h^ = g a Pdl p h^ = F^{h){dh,dh), V / u, l / = 0,...,3, 

where g^ v = + h^ iu and the nonlinear term F is as in Lemma WM. We assume that the initial data 
(g,dtg)\t=o = (5o><?i) are smooth, the Lorentzian metric is of the form 

go = —a 2 dt 2 + goijdx l dx^ 

and 

1) obey the wave coordinate condition 

(2.2) g aa 'd a g a/fl = ^g aa 'd^g aQ/ , V M = 0, ...,3, 

9 Since the initial metric is always of the form gij — (1 + 4M/r)Sij + o(r _1 ) with M > 0, data coinciding with the 
Schwarzschild outside of a compact set is the closest analogue of compactly supported or rapidly decaying data usually 
considered in small data global existence results for nonlinear wave equations. 

10 In what follows we shall introduce the reduced wave operator D s = g d^g and note that in wave coordinates 
□ s = Dg, where O g <f> = \g\ ^ 1 ^ 2 d a [g a/3 \g\ 1 ^ 2 dp(f>) is the geometric wave operator 
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2) satisfy the constraint equations 

Ro - M 2 + (trk ) 2 = 0, V j k 0ij - Vitrko = 0, Vi = 1, ...,3, 

where Rq is the scalar curvature of the metric goij, and the second fundamental form (ko)ij = —l/2a gnj. 

3) we assume that the metric (go)ij coincides with the spatial part of the Schwarzschild metric g s 
(in wave coordinates): 

(9oh = -^TI7 dr2 + ( r + 2M) 2 (d6 2 + sin 2 9 d<p 2 ), r > 1 
r — 2M 

and g\ = for r > 1. Moreover, we assume that the lapse function a 2 (r) = (r — 2M)/(r + 2M) for 
r > 1 and a{r) = 1 for r < 1/2 

4) the data (ho, hi) = (go — m,gi) verify the smallness condition 



(2.3) e = ^E N {0) +M<e , 
where N > 10 and 

(2.4) E N (t) = sup V \\dZ J h(r, IWl, 

Here Z 1 is a product of \I\ vector fields of the form d{, Xidj — xjdi, tdi + X{dt and tdt + x l di. Then 
there exists a unique global smooth solution g with the property that for some constant Cn, 

(2.5) E N (t) < 16e 2 (l + t) 2CNE , 

\\g^(t) - m^|| L » < C N e{l + t)- 1+CN£ . 

Remark 2.2. The existence of data satisfying the assumptions of the theorem follows from the work of 
|Co| . |C-D| . as argued in section 4. 

A corollary of the above result is the global stability of Minkowski space for a restricted set of 
initial data. 



Theorem 2.3. Let (M. , go, kg) be the initial data set for the Einstein vacuum equations R^ v = 0. 
Assume that relative to some system of coordinates (x\, X2, x$) the metric go coincides with the spatial 
part of the Schwarzschild metric gs outside the ball of radius one, 

go = (l + —fdx 2 , r > 1, 
r 

while the second fundamental form ko vanishes for r > 1. In addition, we assume that relative to that 
system of coordinates go, M and ko satisfy the smallness condition 

Yl W dI x(9o- S)\\ L 2 (Bl) + \\ dI xk \\ L 2(B 1 ) + M < e. 

0<\I\<N 0<[J|<JV— 1 

Then there exists a future causally geodesically complete 11 solution g together with a global system of 
wave coordinates with the property that the curvature tensor of g relative to these coordinates decays 
to zero along any future directed causal geodesic. 

11 For the definition see H-E and section 16 of this paper. 
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We now outline the strategy of the proof. 

Remark 2.4. Throughout the paper we shall use the notation A < B for the inequality A < CB with 
some large universal constant C. In our estimates we will make no distinction between the tensors 
h a p = g a p — m a p and H a p = rn aa imppi{g a ^ — m a P), since H = —h + 0{h 2 ) and the terms quadratic 
in h are lower order. 

The continuity argument For the proof we let 5 be any fixed number < 5 < 1/2. Let g be a 
local smooth solution of the reduced Einstein equations (|2.1|) . We start with the weak estimate 

(2.6) E N {t) <Ue 2 {l + t) 25 

By assumptions of the Theorem the estimate Q2.6JI holds for t = 0. Let [0, T] be the largest time 
interval on which (|2.6fl still holds. We shall show that if e > is sufficiently small then on the interval 
[0, T] the inequality ()2.6|) implies the same inequality with the constant 64 replaced by 16. It will then 
follow that the solution and the energy estimate (|2,6|) can be extended to a larger time interval [0, T'] 
thus contradicting the maximality of T. This will imply that T = oo and the solution is global. We 
will in fact prove that for a sufficiently small e the stronger estimate 1)2.5(1 holds true on the interval 
[0,T]. 

The global Sobolev inequality of Proposition 19.21 and the weak energy estimate ()2.6Jl imply the 
pointwise decay estimates: 

,2.7) E l»*'»ft.)l< {1 + t + Z + + t r|)V» - ^ 

|I|<jV-2 V U 

From the assumption that the constant 5 < 1/2 we derive the following weak decay estimates 

(2.8) {dZ^i^x)] <C7e(l + t + r)- 1 / 2 - 7 (l + |t-r|)- 1 / 2 - 7 , \I\ < N - 2 

with some fixed constant 7 > 0. The weak decay estimates (|2.8|) will lead to much stronger decay 
estimates in Theorem 114.11 In turn, using the stronger decay estimates in Theorem 1 1 4 . 1 1 we will be 
able to obtain stronger energy estimates in Theorem ll5.11 i.e. ()2.5|) . These in particular will enable us 
to show that the estimate (j2.6|) holds globally in time and conclude the proof. We remark that in the 
course of the proof all constants will be independent of e > but they will depend on a lower bound 
for 7 > (and hence on an upper bound for 5 < 1/2). 

As described above, the proof is a direct consequence of three results. First is the global Sobolev 
inequality of Proposition 19.21 introduced by S. Klainerman |Klj . giving decay estimates in terms of 
energy estimates for the generators of the Lorentz group. The second ingredient is the improved decay 
estimates in Theorem 114.11 The final component is the energy estimates in Theorem 115.11 which rely 
on the improved decay estimates. 

Weak decay estimates. As pointed out above we may start by assuming the weak decay estimate 
(j2,8() . Furthermore, since the solution g = m + h coincides with the Schwarzschild solution of mass 
M < e in the region r > t + 1, we have 

(2.9) |Z 7 /i(i,x)[ < e(l + r + ty 1 , when |x|=t + l 
Hence integrating (|2.8|) from the light cone, where (|2.9|) holds, we get 

(2.10) |Z 7 /i(i, x)\ < e(l + r + t)- l l 2 ^{l + \t- r]) 1 / 2 ^, 
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Since the vector fields span the tangent space of the outgoing light cones r — t = q we infer, with d 
denoting the derivatives tangential to the cones, that 

(2.11) \BZ I h\ < e(l + r + t)- 3 / 2 -^(l + \t- r\) l l 2 -\ 

This means that, close to the light cone t = r, derivatives tangential to the forward light cones decay 
quite a bit better than the expected decay rate from 1)2. 8j) for a generic derivative. 

Wave coordinate condition. As we shall see below certain components of the tensor h decay 
faster than others. This can be seen upon introduction of a null frame of vector fields L = dt + d r , 
L = d r — dt and Si,S2- two orthonormal vectors tangential to the sphere of radius r in R 3 . The 
first improved estimates come from the wave coordinate condition 1)2. 2j) . Writing g a p = m a /3 + h a p we 
obtain from ()2.2|) that 

m^dji^ = d^m a ^h a p + 0(hdh) 

In particular, contracting with a vector field T € T = {L, Si, S2} and using that for any symmetric 
2-tensor k, m al3 k a {3 = —kiL + § kAB, implies that we can express the transversal derivative 3l of 
certain components of h in terms of the tangential derivatives that decay better and a quadratic term 

\{dh) LT \ < \dh\ + \h\ \dh\ <e(l + t + r)- 1 - 2 ^, \h LT \ < e(l + \t - r|)(l + t + r)" 1 

Even though the estimate above does not give a better decay rate for all components of h it gives 
the decay exactly for those components which, as it turns out, control the geometry, i.e., they lead to 
stronger energy and decay estimates. 

The above estimates will be sufficient to obtain improved estimates for the lowest order energy of 
h. However, in order to get estimates for the energy of Z^h we commute the vector fields Z through 
the equation for h. This generates additional commutator terms. The main commutator terms are 
controlled with the help of the following additional estimate from the wave coordinate condition: 

(2.12) \{dh) LT \ + \(dZh) LL \ <e(l + t + r)" 1 "^, \h LT \ + \{Zh) LL \ < e(l + \t- r|)(l + t + r)' 1 

We now describe derivation of the stronger decay and energy estimates. 

Stronger decay estimates. We rely on the following decay estimate for the wave equation on a 
curved background 12 : 

(2.13) ||(1 +t + r)80(t, -)\\ L oo < C f (1 + r)\\U g ^{r, dr 

Jo 

+ C sup ^||Z 7 0(r,-)||L- +C [ ^(l + rJ-^IZ^Cr.OIUcodr 
o<-r<t mi Jo |7| < 2 

The estimate ()2.13)l will be applied to the components of the tensor h. The term Z 1 h on the right 
hand side of the estimate will be controlled with the help of the weak decay estimates, and thus the 
decay rate of h will be determined in terms of decay of U g h. The L 00 — L°° estimate 1)2.13)1 does not 
rely on the fundamental solution as does the more standard L 1 — L°° type estimate. This estimate 
was used |Llj in the constant coefficient case and here we establish it in the variable coefficient case 

12 Recall that the reduced wave operator IZ1 8 = g al3 d\jj. 
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only under the assumption of the weak decay of all of the components of the metric g and the stronger 
decay of the components of g controlled by the wave coordinate condition. This analysis is by itself 
very interesting but we will not go into it here and just refer the reader to the following sections. 

We now analyze the inhomogeneous term in the equation for h^. The tensor hav = g^ u — 
verifies the reduced Einstein equations of the form: 

Ugh^ = F IJ , 1/ (h)(dh,dh), 

(2.14) F^(h)(dh,dh) = P(d^h,d v h) + Q^(dh,dh) + G^{h){dh,dh), 

(2.15) P(d, t h, d v h) = ^tThd v txh - \d^d v h a ^ 

Here are linear combinations of the standard null- forms and G^ u (h)(dh,dh) is a quadratic form 
in dh with coefficients as smooth functions of h vanishing at h = 0. The weak decay estimates imply 
that the last two terms decay fast 

(2.16) \Q^(dh, dh)\ + \Gp,(h){dh, dh)\ < \dh\ \8h\ + \h\\dh\ 2 <e 2 {l + r + i)- 2 ~ 2 T(l + \t - r\)~^ 

The problematic term is P(d^h, d v K) since a priori the weak decay estimates only give the decay rate 
of e 2 (l + r + t) _1_27 (l + \t — r|) -1-27 , which is not sufficient in the wave zone t ~ r. The crucial 
improvement comes as a result of a decomposition of the tensor P(d^h,d u h) with respect to a null 
frame {L, L, S±, 82}- Let T G T = {L, £1, £2} be any of the vectors generating the tangent space to 
the forward Minkowski light cones and U E U = {L, L, S\, S2} denote any of the null frame vectors. 
Define, for an arbitrary symmetric two tensor k, \k\ru = J2reT ueu = iT^U^k^ul. It then follows that 



(2.17) \P(dh,dh)\ TU < \dh\\dh\<e 2 (l + r + t)- 2 ^(l + \t-r\)-'^ 
On the other hand, the absolute value of the tensor P(dh, dh) obeys the estimate 

(2.18) \P(dh, dh)\ < \dh\ 2 TU + \dh\ LL \dh\ 

We now decompose the system of equations for h with respect to the null-frame 

(2.19) \U g h\ru <e 2 (l + r + tr 2 ~ 2 \l + \t- r\)~ 2 \ 

(2.20) \D g h\uu < \dh\ 2 ru + e 2 (l + r + t)" 2 " 2 ^! + \t - r\)- 2 \ 



where in the last inequality we also used the improved decay estimate for dhu obtained from the 
wave coordinate condition. The result is a system of equations where the components n\ g IiTu have 
very good decay properties, while \D g huu f° r the remaining non tangential component depends, to the 
highest order, only on the components Htu satisfying a better equation. An additional subtlety in the 
above analysis is the fact that contraction with the null frame does not commute with \D g (or even 
with □). However, the decay estimate (|2.13|) for the wave equation only uses the principal radial part 
of □: d 2 — r~ 2 d 2 — 2 r~ 1 d r , which respects the null frame. This analysis results in the improved decay 
estimates 

(2.21) \dh\ TU < Ce{\ + t)-\ \dh\ < Ce{\ + t)^ 1 ln(2 + t) 
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The energy estimates. We rely on the following energy estimate for the wave equation, which 
holds under the assumption that the above decay estimates hold for the background metric g: for any 

7 > 

,2 , r f 7 I^P _ f , nMa 



(2 - 22) L m+ L h A^\t-r\)^ ^k m+C£ h J,— H \L LP^ 

This implies that the energy of a solution of the homogeneous wave equation O g <p = grows but at 
the rate of at most (1 + 1) . The presence of an additional space-time integral containing tangential 
derivatives on the right and side of H2.22|) is crucial for our analysis. This type of estimate in the 
constant coefficient case basically follows by averaging of the energy estimates on light cones used e.g. 
in |S1| . We also note that the energy estimates with space-time quantities involving special derivatives 
of a solution were also considered and used in the work of Alinhac, see e.g. |A2j . |A3j ) . In our work 
we use the space-time integral with derivatives spanning the tangent space to outgoing light cones and 
weights dependent on the distance to the Minkowski light cone r = t + 1. We emphasize that the 
energy estimate (|2.22|) is proved only under the assumption of the weak decay of all components of 
the background metric g together with the strong decay of the components controlled from the wave 
coordinate condition. 

It is worth noting that a combination of the energy estimates of the type (|2.22j) and Klainerman- 
Sobolev inequalities would also yield a very simple proof of the small data global existence result for 
semilinear equations = Q(d(j),d(p) obeying the standard null condition. This fact appears to be 
previously unknown. 

The energy estimate ()2.22(l will be applied simultaneously to all components of the tensor h. As in 
the equations (|2.19|) . (|2.20|) the inhomogeneous term obeys the following estimate: 

\O g h\ <e(l + r + t)- 3 / 2 ^(l + \t- r\f/ 2 ^\dh\ + e{\ + ty^dh], 

where in the last inequality we used the improved decay estimate for the \dh\-ru components. The 
energy estimate (jlOj) will then imply that 



E Q (t) < 16e 2 (l + t 



Higher order energy estimates. In addition to the energy estimates for the components of the 
tensor h we need estimates for the higher vector field derivatives of h: Z 1 h with Minkowski vector 
fields Z = {d a , S}. To obtain these estimates we apply Z 1 to the equation D g h^ u = for h. 
Applying vector fields to the nonlinear terms F^ u yields similar nonlinear terms for higher derivatives 
and these are can be dealt with using the estimates already described above. We must note however 
that this is where the additional space-time integral involving the tangential derivatives on the left 
hand side of the energy estimate ()11.3|) becomes crucial. Consider for example the term dh ■ dZ l h 
generated by one of the null forms in F^y. We estimate its contribution, with the help of the weak 
decay estimates, to the energy estimate as follows: 

m is* < < + Cem ' h? 



(1 + i)l/2+ 7 (1 + | t _ r |)l/2+7 - (1+4)1+27 (1 + \ t - r |)l+2 7 

The integral of the first term is easily controlled by the energy on time slices times an integrable factor 
in time. The space time integral of the second term is in fact part of the energy (|2.22|) . and if we 
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choose e sufficiently small this term can be absorbed by the space time integral on the left. The idea 
with the space-time integral is that one can use the extra decay in \t — r\ when one does not have full 
decay in t. 

The more serious problem in higher order energy estimates lies however in the commutators between 
Z 1 and the principal part O g = g al3 d a dp. 

The commutators. Writing g a/3 = m a ^ + H a ^ with H a/3 = —m aa mr" h a ipi + 0(h 2 ), we show 
the following commutator estimate 13 

(2,3) |[z,q,M < + lZ Z]:lT ) £ T^T-r E V*'« 

by the weak decay assumptions (|2.10|) and the improved decay from the wave coordinate condition 
(|2.12|) . We should note that for a generic quasilinear wave equation commutators with Minkowski 
vector fields Z give rise to uncontrollable error terms. In the special case of the equation U<p = 
this problem can be overcome by modifying the vector fields Z, |A3j . In our case it is the wave 
coordinate gauge that provides additional cancellations. 

This commutator estimate applied to <fi = h a s together with the analysis in the previous section now 
gives estimates for the energy E\ as well as for the stronger decay estimates for the second derivatives 
of h, Q2.26JI with | J\ = 1. This commutator will also show up as a top order term [Z, O g ] ■ Z I ^ 1 h a p in 
the energy estimate for Z 1 h and the resulting term can be dealt with in the same way. 

The other top order term generated by the commutators [Z I ,D g ]cf) is of the form (Z 1 H a ^)d a dp. 
We first apply the pointwise estimate 

\Z I H\ LL ^ ^ lar7K 

\K\<1 



l(^«,M<<;( i LJ- + J T ij £ \az 



To deal with its contribution to the energy estimate we use the Poincare estimate with a boundary 
term 

f ^Hlfrdx f , r f \d r Z T H\ 2 r r dx 

(2 ' 24) J (l + \t-r\)^ - C J {Z HlLLdS + C J [l + \t-%** ' a> " 1/2 ' ° + 

R 3 ' S (t+1) i?3 

together with the fact that h is Schwarzschild outside the cone r = t + 1, because of the inward bending 
of the Schwarzschild light cones, and hence there \Z I h\ < Ce/(1 + t). The way coordinate condition 
implies that \dZ I H\j J ] J can be estimated by \dZ 1 K\ and lower order terms. The term involving \dZ I H\ 
is then controlled by the space-time integral on the left hand side. 

One can use a similar but more trivial argument for decay estimates, i.e. 

\Z T H\ LL < |Z J # LL | r=m + (1 + \t - t^Z 1 H LL \ LOC 

The lower order terms. So far we have only discussed the top order terms, but there will also 
be several lower order terms (relative to |7| = k + 1) to deal with. These are typically of the form 

(2.25) \dZ J h\ \8Z K h\ or \Z J h\ \d 2 Z K ~ l h\ < \dZ K h\ 



13 This commutator estimate applies to the vector fields Z — {d a , Qap}- For the scaling vector field Z — S = x a d a the 
commutator expression should have the form D g S — (S + 2)D g . 
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with | J\, \K\ < 1 1\ = k + 1 The lower order terms are dealt with using induction. We describe the 
induction argument for the decay estimates. From this it will be clear how it also proceeds for the 
energy estimates. We will inductively assume that we have the bounds: 

(2.26) \dZ J h\ + \Z J h\{l + \t-r\)- 1 <C k t- 1+Ck£ , \J\ < k 

The terms in (|2.25j) can then be estimated by C 2 e 2 t~ 2+2Ck£ . Including the top order terms using (|2.23[) 
applied to <fi = Z I ~ 1 h, and using (|2.13|) applied to t} g Z I h we get an inequality of the form 

* C eM(s) Ce 2 
l + s + (l + s) ] 



( 2 - 27 ) M (t) < I Ai-ce ds 



where M{t) = (1 + The Gronwall's inequality then gives the bound M(t) < C(l + 

t) 2C£ . 

3 The Einstein equations in wave coordinates 

For a Lorentzian metric g^u, where fj,, v = 0, 3 we denote 

(3.1) T* v = ^g xs (d^gsu + d v g 5lJi - d s g^) , 

the Christoffel symbols of g and 

( 3 - 2 ) R^vS = fyF^ - d v T^ s + T p gT^v - Y* V Y£ 5 

its Riemann curvature tensor with R^ v = R fl a ua , the Ricci tensor. 
We consider the metric g satisfying the Einstein vacuum equations 

(3.3) = 0. 
We impose the wave coordinate condition: 

(3.4) r A :=^r a ^ = o 

It follows that assuming (|3.4[) we have that the reduced wave operator \3 g = g a/3 

(3.5) n g = U g = -J^dag^y^dp 

The following lemma provides the description of the Einstein vacuum equations in wave coordinates 
as a system of quasilinear wave equations for g^ u . 

Lemma 3.1. Let metric g satisfy the Einstein vacuum equations (|3.3|) together with the wave coordinate 
condition (|3,4[) . Then g„v solves the following system of reduced Einstein equations: 

(3.6) Ogg^ = Pidf.g, d u g) + Q^{dg, dg) 
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where 

(3.7) P(d^g, d v g) = \g aa '^<w - -^V^-W d v g a ,p, 

(3.8) Q^u{dg, dg) = d a g^ g aa ' g w ' d a >gpi v - g aa g pp (d a gp^ dp>g a / u - dpigp^ d a g a >„) 

+ 9 aa ' 9 m ' {d^g a >p>d a gp v - d a g a > p>d^gp u ) + g aa ' g W (d u g a >p>d a gp fl - d a g a >p>d u gp fl ) 

+ ^9 aa ' 9 W ' (d/3'9aa'dfj,gp u - d^g^d^gpv) + -g aa ' g m ' (d^g^d^gp^ - d v g aa! d^g^ 
Furthermore, the wave coordinate condition (|3.4[) reads 

(3.9) g aP d a g^ = -g^d^p, or d a g av = -g^d^ 

Proof. The proof of (j3,9j) is immediate. 
We now observe that 

dag /3 v = + T afJi/ 3, where r Mai/ = g a \T* v . 

It follows that g a \dpT^ u = dpT mv - (Tp a>l + Yp\ a )T* v so also usin S that r «A/3 = Fpxa we obtain 

(3.10) R^iavP = 9a\R^ u f3 = dpT — d u T ^p + T v \ a T ^p — TaxpT^ v 
It follows from ((221) that 

(3.11) g al3 (d^d a gp v - ^d^d v g a p^ = -d^g al3 (d a gp u - ^d u g a p^j = g aa> ' g^' d^g al p,{d a gp v - ^d u g a p 

and hence 
(3.12) 

gO/3 

g ap (daT^pv - dyT^pa) = — {d a d ll gp v + d a d u gp^ - d a dpg^) — (d^gpa + d v d a gp^ - d v dpg m ) 

= — ^-dadpg^ + —(dad^gpu + d v dpg m - d v d^gp a ) 
= -^g^dadpg^ + ^g aa ' g^' {d^g a > p>d a gp u + d v g a > p<d a gp^ - d v g a > p'O^p) . 
Here by (|3.9|) we can write 

(3.13) g aa ' 'gPP' d^g^ p<d a gp u = g aa g^ d a g a < p>d^gp u + g aa g p/3 {d^g a > p>d a gp u - d a g a > p'd^gp u ) 
= 7i9 aa 'g PI3 ' d/3>ga'a d^gpy + g aa ' ' g m ' [d^p' d a gp u - d a g a >p> d^gp v ) 



g aa' g pp> Q^g a , a dp>gp u +g aa ' g P13 ' {-(dp>g a i a d^gpy-d^g^a dp l gp v ) + {d IJ ,g a >pi d a gp u -d a g a > p> d^gpy 



2 

aa' g i3{3' Q^g a , a d v gpp<+g aa ' g 1313 ' (~(dprg a > a d^g^-d^a dp'gp u ) + (d fl g a/ p / d a gp u -d a g a > pi d^gp v ) 
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Hence by ()3.13j) and (|3.13|) with fi and v interchanged we get 



(3.14) 

aa' g pp' [Q^p, d a gp v +d v ga>pi d a g^-d v g^p< d^g a ^j = g aa ' g pp ' {-d^a d v gpp>-- 

1 — „' aa> 



+ -j9 aa ' 9 PP ' ((pn9a'p' d a g/3u ~ d a g a 'p> d^gp u ) + {d v g a >p d a g^ - d a g a >p> d u gp^j 

-g aa ' gW (^(dpga'a d^gp u - df,g a > a dp>gp u ) + (dp>g a > a d v g^ - d v g a > a dp/gp^ 

On the other hand 

(3.15) TuapT™ 13 = - (d u gp a + dpg au - d a gp v )g aa ' ' g w ' (c^'a' + dp>g a ^ - d^gp^) 

= jdvQaP g aa ' g w ' d^p, + -0 a 0/9/i 9 aa ' 9 PP ' 9 a '9p'u ~ i^d^gp^ g aa ' g PP ' Qp'9oj v 
= 9 aa 9^ {jdv9 a p df.ga'p' + -jd a gp^ d a 'gp> u ~ -^P'9Pii 
- ^9 aa ' 9 PP ' (p a 9Pn dp,g a i u - dpg^ d a g a i u ^j 

= 9^ ' g PP ' (^dv9aP 9^9a'P' - g^5/3/3' 9v9aa> + ~^pa9Pn d a '9p'i 



- ^9 aa 'g^' (pagpn dprg a > v - dprgp^ d a g a > u 

where the last inequality follows from ()3.9|) . 

Taking the trace of (|3.1()|) and using (|3.12j) . (|3.4|) we obtain 

(3.16) R^y = -]^g aP d a dpg^ v + T va pT ^ + ^g aa ' gW [d^p, d a gp v + d v g a < p, d a gp^ - d u g a ,p, d^g a p 
Using (|3.15j) and 1)3. 1 4 j) we get 

(3.17) = --g^dadpg^ + g aa ' g^' ( - -d„g a p d ll g a >p> + -d^gpp* d v g aa! 

+ -j9 aa ' 9^ da9Pn d a ,gpi v - ^g aa ' gW {d a 9py, dp>g a > u - dpigp^ d a g all 



1 — „' aa> 
— ( 

T 

1 ™' ani . 



+ T.g aa ' g^' ( ip^ga'P' d a gp v - d a g af p> d^gp u ) + (d v g a >p> dagpf, - d a g a >p' d v g 0ll 



_ g aa' g pp' ^d /3 ,g a , a d^gp u - d^g a , a dp,gp v ) + (dp>g a i a d v gp^ - d u g a , a dp,gp^j 

The result now follows. □ 
Let m denote the standard Minkowski metric 

moo = -1, mu = l, if i = l,...,3, and = 0, if n^v, 

Define a 2-tensor h from the decomposition 
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Let m^ v be the inverse of m^ v . Then for small h 

H nv = g nv _ m nu = _ h m> + o^(h 2 ), where F = m^'m^'^v 

and O^ih?) vanishes to second order at h = 0. 
As a consequence of Lemma 13. II we get: 

Lemma 3.2. If Einstein's equation's (j3.3j) and the wave coordinate condition 1)3. 4 Jl hold then 

(3.18) n g h tMU = F t ,v(h){dh,dh) 

where F^ v {h){dh,dh) is a quadratic form in dh with coefficients that are smooth functions of h. More 
precisely, 

(3.19) F^(h)(dh,dh) = Pid^dvh) + Q„ v {dh,dh) + G^(h)(dh,dh) 
where 

(3.20) P^h, 8 u h) = \m aa ' d^h aa , m^'d v hppi - l -m aa ' m^' 8^ h a p d v h a ipi 
and 

Q^dh^dh) = d a h/3[t m aa ' mP 13 d a ,hpi v - m aa 'm^ ' (d a h^ dp>h a , v - dpihp^ d a h a , u ) 

+ m a01 mP P ' [d phot p d a hp v - d a h a ipi d^hp v ) + m aa ' ' mP p ' [d v h a ipi d a gh PvL - d a h a ipi d u h Pp ) 

+ -m aa (dpih aa i d^hp u - d^h aa i dp<hp v ) + -m aa ' mP®' {dp>h a0l i d u hp^ - d u h aa , dpihp^ 

is a null form and G u U (h)(dh, dh) is a quadratic form in dh with coefficients smoothly dependent on h 
and vanishing when h vanishes: G flu (0)(dh, dh) = 0. 
Furthermore 

(3.21) m^d a hp^ = ^m^d^h a p + G^h)(dh), or d a H av = ^(m'" + W^d^H * 

where G tl {h){dh) is a linear function of dh with coefficients that are smooth functions of h and that 
vanishes when h vanishes: (7^(0) (<9/i) = 0. 

Observe that the terms in (|3.2U|) do not satisfy the classical null condition. However the trace 
m^h^ satisfies a nonlinear wave equation with semilinear terms obeying the the null condition: 

g^dadpm^h^ = Q(dh,dh) + G{h){dh,dh). 
4 The initial data 

In this section we discuss the initial data for which the results of our paper apply. We shall consider 
the asymptotically flat data, satisfying a global smallness condition, with the property that it coincides 
with the Schwarzschild data outside the ball of radius one. 
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We start by showing the existence of such data. Let (go, ko) be asymptotically flat initial data for 
the Einstein equations consisting of the Riemannian metric go and a second fundamental form ko- The 
initial data for the vacuum Einstein satisfy the constraint equations 

(4.1) i? -(tr£;o) 2 + |/co| 2 = 0, 

(4.2) V j k 0ij - V.trfco = 

We restrict our attention to the time-symmetric case Ro = ko = 0. Then, if (go, ko) is sufficiently close 
to the Minkowski data and go satisfies the parity condition go(x) = go(—x), by the results of Corvino 
|Coj and Chrusciel-Delay |(J-D| one can construct a new set of initial data (g, k) with the properties 
that 

• The initial data (g, k) coincides with (go, ko) on the ball of radius 1/2. 

• (g, k) is exactly the Schwarzschild data (g$,0) of mass M outside B\, the ball of radius one. 
At this point we specify the smallness conditions: 

(4-3) M<e, Yl (\\ d i(9-S)\\ L 2 {Bl) + £ \\d J x k\\ L 2 [Bl) ) < e 

0<|/|<7V 0<| J"|<A^— 1 

for some sufficiently large integer N. Here d], denotes the derivative d^.-.d^, where (I\, .., I n ) is an 
arbitrary multi- index with the property that 1\ + .. + I n = \I\. 

We have two expressions for the Schwarzschild metric in isotropic and wave coordinates: 

(4.5) g s = - r -^Ldt 2 + ^±^-dr 2 + (r + 2M) 2 (d0 2 + sin 2 0# 2 ) 
V ; y r + 2M r - 2M V ' V Y ' 

The expressions g$ and g^ will denote the spatial parts the Schwarzschild metric in respective coordi- 
nates. Observe that 

(4.6) g s = m + —(dt 2 + dx 2 )+0(r- 2 ), 

r 

We now find the coordinate change transforming the metric gs into g s . Set 

M 2 

(4.7) t = r, r = p+ 

P 

In the coordinates r, p the metric g s takes the form gs- This change of coordinates is one-to-one for 
the values p > M. Since the mass M « 1 we can define the change of coordinates r = &(p), where 
$ coincides with the map (|4.7|) for p > 1 and the identity transformation for p < 1/2. Thus we have 
constructed the initial data (g, k) such that 

• The initial data (g, k) coincides (in new coordinates) with (go, ko) on the ball of radius 1/2. 

• (g, k) is exactly the Schwarzschild data (g*,$) outside the ball of radius one. 
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• Moreover, the new data still obeys the smallness condition (|4.3[) . 



The constructed metric is already in wave coordinates on its Schwarzschild part. We now describe 
the procedure which produces the initial data (g, dtg) associated with (g, k) and satisfying the wave 
coordinate condition. 

Recall that a priori we are only given the spatial part of the metric gij together with a second 
fundamental form kij. We now define the full space-time metric g a p on the Cauchy hypersurface £o 
as follows: 

(4.8) g 0i = 0, g 00 = -a(r), 
where the function 

, . r-2M 
a(r) = 7T2M' ^ r > 1 

a(r) = 1, for r < — 

Thus defined metric coincides with the full Schwarzschild metric g s for r > 1. We further define 

(4.9) d t gij = -2ahj 

It remains to determine dtgoa- We find it by satisfying the wave coordinate condition 

g^d^ = -g^dag^ 

Setting a = we obtain 

^g 00 d t goo = -g^digop + -g ij d t gij 
This defines dtgoo- On the other hand setting a = i we obtain 

g 00 d t goi = -g pj d,g i(3 + \g^d l9 ^ 

This determines dtgoi- Observe that since the metric g coincides with the Schwarzschild metric g s , 
already satisfying the wave coordinate condition, outside the ball of radius one, we have that on that 
set the initial data takes the form (g s ,0). Hence we constructed the initial data (g,dtg) with the 
properties that 

• The initial data (g, dtg) corresponds to the initial data (g, k) prescribed originally. 

• (g,d t g) is exactly the Schwarzschild data (g s ,0) outside the ball of radius one. 

• The initial data verifies the wave coordinate condition. 

• The initial data satisfies the smallness condition 

(4-10) Yl (\\ d i(9-m)\\ L 2 {Bl) + Yl R J ^ll Wl) ) < e 

0<|/|<7V 0<\J\<N-l 
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Now with the initial data (g,dtg) we solve the reduced Einstein equations ()3.6|) , It follows from the 
proof of Lemma 13. II that, in the notation T x = g a ^^^, the reduced Einstein equations can be written 
in the form: 

(4.11) R a/3 - \{D a T, 3 + D 3 T a ) - T a N^(g, dg) = 0. 

Here D denotes a covariant derivative with respect to the space-time metric g and N°g are some given 
functions depending on g and dg. Observe that the initial data (g, dtg) were chosen in such a way that 
the wave coordinate condition T x = is satisfied on the initial hypersurface T,q. We now argue that 
this condition is propagated, i.e, the solution of the reduced Einstein equations (|4.11j) obeys T A = 
on any hypersuface Sj. We would have thus shown that a solution of the reduced Einstein equations 
is, in fact, a solution of the vacuum Einstein equations. 

To prove that T A = we differentiate (|4.11jl and use the contracted Bianchi identity D^~R. a p = 
^Z) Q R 

= 2(Z^R Q/3 - ~£> Q R) = D^D a T + D^D p V a - D a D^ - 2D^(T a N^) - D a (T a Nf) 
= D^DpTa + R a7 r 7 - 2{D p Y a )Nf - {D a T a )Nf - 2T a (D (3 Nf) - T a (D a Nf) 

Therefore, T x satisfies a covariant wave equation, on the background determined by the constructed 
metric g, with the initial condition T x = 0. It remains to show that D t T x = on So and the conclusion 
that r A = will follows by the uniqueness result for wave equation. 

We recall that the initial data (g, k) verifies the constraint equations (|4.1[) . (|4.2I) . which imply that 
on Eo 

R TT + ^R = 0, Rt; = 0, 

where T = — (goo) 1 ^ is the un it future oriented normal to Therefore returning to (|4.11|) we 
obtain that 

= R 00 + ^R = -(goo)' 1 D t T + D'Ti, 
= Koi = ^D t Ti + -AT 

This finishes the proof that T x = 0. 

We also know that the time-independent Schwarzschild metric g s is a solution of the Einstein 
vacuum equation R Q/ g = 0. Moreover, since g s satisfies the wave coordinate condition it also verifies 
the reduced Einstein equations (|4.11|) . Since the initial data (g, dtg) = (g s , 0) outside the ball of radius 
two, constructed solution will coincide with the Schwarzschild solution in the exterior of the null cone 
developed from the sphere of radius one in So- 

We end the discussion of the initial data by comparing the light cones of Minkowski and Schwarzschild 
spaces in the wave coordinates of the Schwarzschild space. 

Lemma 4.1. For an arbitrary R > 2M the forward null cone of the metric g s , intersecting the time 
slice t = along the sphere of radius R, is contained in the interior of the Minkowski cone t — r = R. 
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Proof. The null cone intersecting the time slice t = along the sphere of radius R can be realized as 
the level hypersurface u = of the optical function u solving the eikonal equation 

gfd a u dpu = 

with the initial condition that u = on the sphere of radius R at time t = 0. Because of the spherical 
symmetry of the Schwarzschild metric g s and the initial condition we look for a spherically symmetric 
solution u = u(t,r). The eikonal equation then reads 

r + 2M 2 r-2M 2 

— (O t u) = TTWrU) 

r-2M y 1 r + 2M v 1 

Let t = 7(r) be a null geodesic, originating from some point on the sphere of radius R at t = 0, such 
that u(j(r),r) = 0. Then 

dtuj(r) + d r u = 

Substituting this into the eikonal equation we obtain that 

■r + 2M\2 
Jr - 2Ms 

Taking the square root and integrating we obtain that 

r-2M \ 
R-2M, 

Thus the null geodesies are described the curves 

r — 2M 

t = ±[r - R + AM In 



we obtain that 

( r — 2M \ 
r-R + AMln^-^) 



( r — 2M \ 

(r-R + 4Mln ) 

V R-2M) 



R-2M, 

In particular, the forward null cone is contained in the interior of the set t > r — R. □ 



5 The null-frame and null-forms 

Below we introduce a standard Minkowski null-frame used throughout the paper. At each point (t, x) 
we introduce a pair of null vectors (L,L) 

L° = l, £,* = arVI^I. * = 1,2,3, and L° = 1, V = —x i /\x\, i = 1,2,3, 

Adding two orthonormal tangent to the sphere S 2 vectors S\ , S% which are orthogonal to uo defines a 
null frame (L,L, S±, ^2). 

Remark 5.1. Since S 2 does not admit a global orthonormal frame S±, S2 we could alternatively intro- 
duce a global frame induced by the projections of the coordinate vector fields e^. 

Let P be the orthogonal projection of a vector field in R 3 along lo = x/\x\ onto the tangent space 
of the sphere; PV = V — (V,uj}uj. For i = l,2, 3 denote the projection of <9j by 

(5.1) ^~ i^'a,- = Qi - WiO^'a,-, where A\ = (P ei ) j = Sf - uJiJ, i = 1, 2, 3, 

where e, is the usual orthonormal basis in R 3 , and the sums are over j = 1,2,3 only. Let do = L a d a 
and di = $ { , for i = 1, 2, 3. Then a linear combination of the derivatives {Bq, ...,193} spans the tangent 
space of the forward light cone. 
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In what follows A, B will denote any of the vectorfields Si, £2- We will use the summation conven- 
tions 

X A A a = Xf>S lp S? + XPS 2P S%, X A Y A = X a Y?S la S w + X a Y?S 2a S w . 

Obvious generalizations of the above conventions will be used for higher order tensors. 

We record the following null frame decomposition of a vector field X = X a d a : X a = X L L a + 
X—L a + X A A a . Relative to a null frame the Minkowski metric m has the following form 

m LL = rriLL = m LA = uila = 0, m L L = tull = -2, m AB = S AB , 

i.e. m a pX a Y' 3 = -2{X L Y^ + X^Y L ) + X A Y A . Recall that we raise and lower indices of any tensor 
relative to the Minkowski metric m, i.e., X a = m a pXP. We define X Y = m aP X a Y^ = X Q Y a . Then 
Xy = X l Yl + X—Yl + X A Y A . It is useful to remember the following rule: 

XL v" v~L ^ v vA v" 

= — 2 — 2 = X A . 

Then 

m LL = m ^ = m LA = mh A = 0, m L ± = mk L = -1/2, m AB = 5 AB 
i.e. m^X a Y p = -\{X l Yl + XlY l )+X a Y a . 

Definition 5.2. Denote q = r — t and s = t + r the null coordinates of the Minkowski metric m and 
d q = \{d r — dt) and d s = \{pt + d r ), the corresponding null vector fields 



Let k X Y = k a/3 X a YP. Then 

t,r k = m^k^a = - , 
2 



(5.2) tik = m a ^k aP = -l(k L L + kLL)+trk 



where 

(5.3) tr A; = 5 AB k AB = Vk ih and T j = 6 ij - J J 

where the sum is over i,j = 1, 2, 3 only. 
If k and p are symmetric it follows that 

(5.4) p a pk a P = m aa 'mW'p aP k a ,p, 

= j(PLLkLL + VLL^LL + ^VLLkLL) ~ (PAL^BL + P AlM-BL) + S AB S A ' B ' ' VAA'^BB' 

= T {PLhkhL + PLL&LL + ^PLiMll) ~ ^ {PihkjL + Pihkjh) + f^f j Pu'kjj> 

Lemma 5.3. With P(p, k) given by (|3.20|) we have for symmetric 2-tensors p and k: 

(5.5) P(p, k) = ^m^pcpm^kap - \m aa ' m^' p a pk a , p , 

= -^{pLLkLL + PLLk LL ) - -S AB S A ' B ' (2p AA ik B B' ~PABk A ' B > 

+ ^5 AB (2p AL k B L + 2p A L,k BL - p AB k LB - Pll 
i. e. at least one of the factors contains only tangential components. 
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Furthermore 



p aP d a = p Lfi d L + p^d L + p AI3 d A = -\vldL ~ \p[d L + p A/3 d A 



We introduce the following notation. Let T = {L, S\, S2}, U = {L, L, Si, S2}, C = {L} and 
S = {Si, £2}. For any two of these families V and W and an arbitrary two-tensor p we denote 

(5.6) \p\vw= \PPi vPwl l 

(5.7) \d P \vw= £ \{9p)c0rU a VPyn\, 

ueu,V£V,wew, 

(5.8) \dp\vw= £ \{dp) a p~ ( T a V l3 W 1 \ 

TeT,vev,wew, 

Let Q denote a null form, i.e. Q a p{d(j>, dtp) = d a (frdpip — dp<pd a ip if a 7^ f3 and Qo(d(ft,dip) = 
m a Pd a 4>dpip. 

Lemma 5.4. If P is as in Lemma \5. J then 

(5.9) |-P(P,*:)| ^ \p\ru\k\ru + \p\hi\k\ + |p||A:| L l 
If Q(d<f),dip) is a null form then 

(5.10) \Q(d(f>,dn < I^H^I + 19011^1 
Furthermore 

(5.11) |A^d a <W| < (|fc|LL|50| 2 + |fc| I^H^I) 

(5.12) |L a fc a %0| < (\k\ LL \d<P\ + |fc| 1501) 

(5.13) \(d a k al3 )df3(p\ < {\dk\ LL + \dk\) \d<j>\ + \dk\ \dct>\ 

Proof. The proof of (|5.1()j) for the null form Qq follows directly from (|5,2|) . To prove the claim for the 
null forms Q a p use that 

(5.14) di = Li(d s + d q ) + B h i = 1,2,3, d = L (d s -d q ) 

Therefore, 

\Qap(d<l>,di/>)\ = " < C|90| + C|d0| 

The estimates (pTTTJl - tjOSl) follow from (f5~l)) . □ 
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Lemma 5.5. If k a " is a symmetric tensor and 4> a function then 

(5.15) \k aP d a dp4>\ < {\k\ LL \d 2 ct>\ + \k\ \3d<j>\) 
Also, with tr k = b AB k^B = — oJ l oji)kij we have 

(5.16) \k a Pd a d^-k LL dl^-2k LL d s d q 4>-r- l YYkd q 4>\ < \k\ LT \Ba<i>\ + \k\ {\B 2 ^\ + r- 1 \d4>\). 

Proof. The estimate (|5.15|) follow from ()5.4|) . We have 

(5.17) diujj = r _1 ((5ij - uJiUj) = r" 1 ^- 
Furthermore di = di + u>id r , where d r = ui^dj so [<9j,<9 r ] = (BiUJ k )dk and 

(5.18) didj = (di + ujid r )(dj + cojd r ) 

= didj + tJi u k djdk + ujjUJ k did^ + uJiUJjd 2 + (diujj)d r + ujj(diUJ k )dk 

= didj + Ui djd r + ujj Bid r + u>iUJjd 2 + r~ 1 5ijd r — r~ x ujidj 

Furthermore 

(5.19) d di = d t (8i + LUidr) = did t + uJid t d r 
Hence 

(5.20) A^d^ = k 00 d 2 + 2k 0i u i d t d r + k ij ujiujjd 2 r + r -1 tr fc d r 

+ k ij didj - r^k^Uidj + 2k 0j Bjdt + 2k lj LOidjd r 

If we substitute dt = d s — d q , d r = d s + d q and identify 

(5.21) k LL = k 00 -2^*^ + u^, k LL = -k 00 + Pw,^, k^ = k 00 + 2k 0t uJi + 0^ 
and 

(5.22) -fc 0j ' + jfetfc^ = V + fc-V = k L j , k 0j + k ij Uj = -k Q j + kJJ = k^ 
we get 

(5.23) hrPd a dp = k LL d 2 q + 2k LL d s d q + kLipl + r^te fc d 9 

+ k ij didj + r _1 tr fc d s - r' l k ij UiBj + 2k L j djd q + 2k L j djd s 

Finally, we can also write 

(5.24) 2k L j 5jd q = k L j dj(uj k d k - d t ) = k L j cu k djd k - k L j djd t + r~ l k^dj, 

since (djLo k )dk = r~ 1 8j. The inequality (|5.16|) now follows. □ 
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Corollary 5.6. Let (ft be a solution of the reduced wave equation H\ g (f) = F with a metric g such that 
H a @ = g a @ — m a @ satisfies the condition that \H L — \ < i. Then 
(5.25) _ 

Ads " ^& dq " tlH 2g^r~ ) dq{r<t>) + 2^1 ~ A ^ + |F ' LT T ^ + ^ ^ + 1^1 + r "Vl) 
w/iere A w = A = S^didj. 
Proof. Define the new metric 

n a/3 
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The equation g a ^d a dp(f) = F then takes the form 

F 



gr p d a 



which can also be written as 

Ucj) + (g a ? - m a/3 )d a r '' 



-2g L k 

Let k a/3 be the tensor k a @ = (g a/3 — m a ^) Observe that 

k a P = (_2 5 i ^)- 1 (^ + 2m Q V-) = (-2/^)~ 1 (H alS + m af3 (2g L ± + 1)) 

= {-2g LL )~ l {H af} + 2m af3 H LL ) 

Thus, 

(5.26) fc L L = 0, fc LT = (-2g L ±Y l H Lr , trfe = (-2 5 L ^)- 1 (t?F + fl^ 

Moreover, < |JEf|, since <? L — = — I and by the assumptions of the Corollary |JJ i — | < |. 

Now using ()5.16)) of Lemma 15.51 with the condition that Ull = 0, together with the decomposition 

□0 = -dU + A^=-(d t + d r )(d r - d t )r<t> + = -d s 5 n£ + A^. 
r r 

we find that the identity + k al3 d a df3(f> = (— 2g L -)~ l F leads to the inequality 

|4d s <V0 + rfoi^ + tr A:9 g( /)+ (2 ff i ^)- 1 rF| < r|A^| + r\k\ LT \dd<t>\ + \k\ (r |d 2 </>| + |^|) 
Finally, identity 1)5. 26|) and a crude estimate \k\ < yield the desired result. □ 

6 The weak null condition and asymptotic expansion of Einstein's 
equation's in wave coordinates 

Let us now first describe the weak null condition. The results of this section appear in |L-R| . Consider 
the Cauchy problem for a system of nonlinear wave equations in three space dimensions: 

(6.1) - □ m = Fi(u,u',u"), i = l,...,N, u = (iti, ...,u N ), 
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where -□ = -d t 2 + £)J =1 d 2 x3 . We assume that F is a function of u and its derivatives of the form 

(6.2) Fi(u, u', u") = a{ k a(i d a Uj 8 p u k + G^u, u', u"), 

where Gi(u, u', u") vanishes to third order as (u, u' , u" )^0 and af a/3 = unless |a|<|/3|< 2and \/3\>l. 
Here we used the summation convention over repeated indices. We assume that the initial data 

(6.3) u(0,x) =su (x) G C°°, u t (0,x) = eu x {x) G C°° 

is small and decays fast as \x\ — * oo. We are going to determine conditions on the nonlinearity such 
that the equation (|6.1jl is compatible with the asymptotic expansion as |x| — > oo and \x\ ~ t 

(6.4) u(t, x) ~ eU(q, s, u)/\x\, where q=\x\—t, s = eln|x|, u = x/\x\, 

for all sufficiently small e > 0. The linear and some nonlinear wave equations allow for such an 
expansion with U independent of s and the next term decaying like s/\x\ 2 , see {Mil IH2j . Substituting 
(|6,4|) into ()6.1|) and equating powers of order e 2 /\x\ 2 we see that 



(6-5) 2d s d q U i = Al K mn (u)(d^Uj)(^U k ), U\ 

where 



s=0 



(6.6) Al k mn (uj) = ^2 a t a ^ a ^^ where u = {-1,uj) and u a = u ai . . .u ak 

\a\=m,\(3\=n 

In fact, Ou = —E~ l d s d q (ru) + angular derivatives and = Cj^d q + tangential derivatives. 

One can show that ()6.1|) - ()6.3|) ^as a solution as long as elogi is bounded, provided that e > is 
sufficiently small and the solution of (|6.5j) exists up to that time, |.T-K| IH1| IH2[ IL1| IL2j . The only 
exception is the case A 3 i00 ^ 0, which has shorter life span. In cases where the solution of (|6.5f) blows 
up it has been shown that solutions of l|6.1j) - (|6.3|l also break down in some finite time T £ < e c ^ £ , 
|Jll IH11 lAlj . John's example was 

(6.7) Du = ut An 

for which (|6.5[) is the Burger's equation (2d s — U q d q )U q = 0, which is known to blow up. The equation 

(6.8) Uu = u 2 t 

is another example where solutions blow up, for which (|6.5|) is d s ll q = U 2 , that also blows up. 
The null condition of K2 is equivalent to 

jk 
~i mn 



(6.9) A{ k mn (uj)=0 for all (i,j,k,m,n), to G S 2 



The results of [CH, |K2j assert that (fO ]l - (fOjl has global solut ions for all sufficiently small initial data, 
provided that the null condition is satisfies. In this case the asymptotic equation (|6.5f) trivially can 
be solved globally. Moreover, similar to the linear case, its solutions approach a limit as s — > oo and 
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the solutions of ()6.1|) - (|6.3j) decay like solutions of linear equations. A typical example of an equation 
satisfying the null condition is 



(6.10) Du = u\ - IV^I 2 

There is however a more general class of nonlinearities for which solutions of (J6.5|) do not blow up: 
We say that a system (|6.1|) satisfies the weak null condition if the solutions of the corresponding 

asymptotic system (|6.5(l exist for all s and if the solutions together with its derivatives grow at most 

exponentially in s for all initial data decaying sufficiently fast in q. 

Under the weak null condition assumption solutions of (|6.5|) satisfy the equation ()6.1j) up to terms 

of order e 2 /\x\ 3 , but need only decay like e/\x\ 1 ~ c ' £ . An example of the equation satisfying the weak 

null condition is given by 

(6.11) Du = uAu 

In |L2j it was proven that ()6.11j) have small global solutions in the spherically symmetric case and 
recently |A3| established this result without the symmetry assumption. The equation (jfi.llj) appears 
to be similar to ()6.7|) but a closer look shows that the corresponding asymptotic equation: 

(6.12) (2d s - Ud q )U q = 

has global solutions growing exponentially in s, see |L2j . The system 

(6.13) Uu = v%, Uv = 

is another example that satisfies the weak null condition. The equation (|6.13j) appears to resemble 
(|6.8|) . The system however decouples: v satisfies a linear homogeneous equation and given v we have 
a linear inhomogeneous equation for u, and global existence follows. The corresponding asymptotic 
system is 

(6.14) d s d q u = {d g v) 2 , d s d q v = o 

The solution of the second equation in (|6.14f) is independent of s: V q = V q (q,u>) and substituting this 
into the first equation we see that U q (s,q,io) = sV q (q,u) 2 so du only decays like |x| _1 ln|x|. 

We show below that the Einstein vacuum equations in wave coordinates satisfy the weak null 
condition, i.e. that the corresponding asymptotic system (|6.5|) admits global solutions. In fact, each 
of the quadratic nonlinear terms in the Einstein equations is either of the type appearing in (|6.1U|) . 
(l6~TTTl or (ffHHJ) . 

Theorem 6.1. Let h be a symmetric 2-tensor and let 

(6.15) hfiu(t, x) ~ ell^ u (s, q, u)/\x\, where q = \x\ — t, s = eln\x\, u = x/\x\. 

be an asymptotic ansatz. Then the asymptotic system for the the Einstein equations in wave coordinates 
(|3.18|) . obtained by formally equating the terms with the coefficients e 2 \x\~ 2 , takes the following form: 

(6.16) {2d s - U LL d q )d q U^ = L^L u P(8 q U, 8 q U), V^x, v = 0, 3 
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where Ull = m aa ' mP^' U a i piL a Lp and P(d q U,d q U) = \d q trU d q trU — \d q U a p d q U a/3 . The asymptotic 
form of the wave coordinate condition (|3.21[) is 

(6.17) 2d q U Lfl = L^dgtrU, V/x = 0,...,3, 

where Ul^ = rn aa U a >nL a and trU = ra a ^U a ^ The solution of the system (|6.16|) - (|6.17jl exists globally 
and, thus, the Einstein vacuum equations ()3.18|) in wave coordinates satisfies the weak null condition. 
Moreover, the component dgUiL grows at most as s while the remaining components are uniformly 
bounded. 

The asymptotic form (|6.16j) follows by a direct calculation from (|3.18|) . Observe that the null form 
Qfiv(dh, dh) disappears after passage to the asymptotic system. 

Next we note that (|6.17|) is preserved under the flow of (|6.16fl . Contracting (|6.16j) with L^L V we 
obtain 

{2d s - u LL d q )d q u LL = o, 

which can be solved globally. More generally, contracting (|6. 16|) with the vector fields {L,Si,S2} we 
obtain 

(6.18) (2d s - U LL d q )dgU TU = 0, if Te{L,Si,S 2 } and U G {L, L, S u S 2 }, 

which can be solved globally now that Ull has been determined. Note that the components d q UTU are 
constant along the integral curves of the vector field 2d s — ULhdq- The remaining unknown component 
Ull can be determined by contracting the equation (|6.16j) with the vector field L. 

(6.19) (2d s - U LL d q )d q ULL = 4P(d q U, d q U), 

By Lemma lS"3*l the quantity P(d q U, d q U) does not contain the term (d q ULL) 2 ■ Thus, the equation (|6.19|) 
can be solved globally and produces solutions growing exponentially in s. A more precise information 
can be obtained from the asymptotic form of the wave coordinate condition (|6.17j) . For contracting it 
with the null frame {L, S\, S2} we obtain d q ULT = 0, if T £ {L, Si, S2}. Therefore, 

(6.20) p(d q u, d q u) = -±5 AB 5 A ' B ' (2d q u AA , d q u BB , - d q u AB d q u A , B ^ - U AB d q u AB d q u L u 

It follows from (|6,18|) that P is already determined and is, in fact, constant along the characteristics 
of the field 2d s — U LLd q . Therefore, integrating (|6.19|) we infer that d q ULL grows at most like s. 

7 Vector fields and commutators 

Let Z G Z be any of the vector fields 

tt af 3 = -x a dfs + xpd a , S = td t + rd r , d a , 

where xq = —t and x% = x l , for i > 1. Let / = (ti, ...,tfc), where = 1, be an ordered multiindex 
of length |/| = k and let Z 1 = Z L1 ■ ■ ■ Z Lk denote a product of |/| such derivatives. With a slight 
abuse of notation we will also identify the index set with vector fields, so I = Z means the index I 
corresponding to the vector field Z. Furthermore, by a sum over I\ + I2 = I we mean a sum over all 
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possible order preserving partitions of the ordered multiindex / into two ordered multiindices I\ and 
J 2 , i.e. if I = (ii,...,i k ), then h = (i h , t in ) and h = (i in+1 , t ik ), where h,...,ik is any reordering 
of the integers 1, ...,k such that i\ < ... < i n and i n +± < ... < it and ix>— >*fc> With this convention 
Leibniz rule becomes Z T (fg) = J2i 1 +i 2 =i(Z 11 f)(Z l2 g). We denote by 9 the tangential derivatives, i.e., 
5 = {5q, di, 82, 83} and note that the span of the tangential derivatives {8q, d\, 82, 83} coincides with 
the linear span of the vectorfields {8 s ,ds 1 ,ds 2 }- 

Lemma 7.1. We have the following expressions for the coordinate vector fields: 

tS — x r £loi 



(7.1) a, 



(7.2) d r = u l di 



t 2_ r 2 ' 

t^Qoi — rS 



t 2 



-x j Q,jj + t£l 0i - XjS _ _ XjS ^ XjX j fl j Q(h 
{ ' { ~ W^7 2 - ~t 2 -r 2 t(t 2 -r 2 ) t 



In particular, 



(7.4) ft = _( ft + ^) 4 = ^-^ r ; 
Lemma 7.2. For any function f we have the estimate 

(7.5) (1 + t + \q\)\3f\ + (1 + \q\)\8f\ < C \Z T f\, \8f\ < \3f\ + I VI 

|/|=l 

where \8f\ 2 = \8 f\ 2 + \8if\ 2 + |<9 2 /| 2 + |<9 3 /| 2 and <9 = <9 S . Furthermore 

(7.6) |a 2 /| < ^ £ r , 

|/|<2 

^ere j<9 2 /| 2 = E«,/s=o,i,2,3 l^/| 2 - 

Moreover, if k a ^ is a symmetric tensor then 

(7.7) 1**^ < C( TT |^ R + £ 1^*1 

Proof. First we note that if r + 1 < 1 then (|7,5|) holds since the usual derivatives 8 a are included in 
the sum on the right. The inequality for \df\ in (|7.5jl follows directly from (|7.4|) : one has to divide into 
two cases r < t and r > t and use two different expressions depending on the relative size of r and t. 
The inequality for \df\ in (|7.5jl follows from Q7.1JI and the first identity in j|7.H)) . 

If i + r < 1 then Q7.6[) follows from (|7.4|) . since [9^1 < Cr _1 and the sum on the right of l|7.6[) 
contains the usual derivatives. Since |fiy-a;fc| < C and Cl^r = tlijt = 0, for 1 < i,j < 3 it follows, by 
applying 9j = r u^fiy- to the expressions in (|7.4|1 . that 

(7.8) |%/| < Cr- 1 ^ + r)- 1 ]T |Z 7 /|- 

|J|<2 
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Once again we distinguish the cases r < t and r > t and use different expressions for <9j. With the 
notation do = 2d s (|7.8f) holds also for a = 0. Since [<9 s ,di] = it only remains to prove (|7.6|) for d 2 . 
Since StJ = 0, \n 0i uJ j \ < Ctr' 1 , S{t + r) = 2(t + r) and \n oi (t + r)\<C{t + r) ,fTSJ) follows also for 
d 2 

The inequality ()7.7j) follows from Lemma [5.51 ()7.5|) and the commutator identity [Z, dj\ = a\dj. □ 
Lemma 7.3. Suppose O g (f) = F. Then 

r\H 



Proof. By Corollary 15.61 



F LL „ tr# + n , , rF 



2^ * 2(^ir r q w^2g L ± 

< r\A^\ + r|H| LT |^| + (r |<9 2 </>| + \8<f>\ + r" 1 

where A w = 5 lJ didj. Here all the the derivatives can be reexpressed in terms of the vector fields Z and 
d q using [721 yielding the expression (|7.9|) , Note that 

- < E|/| = i < E|/|<! I^VI < E m < 2 l^0l 
i + ~ i + t + ~ (i + |?|)(i + t+ 

□ 

Lemma 7.4. Let Z = Z^dn be any of the vector fields above and let Cq be defined by 

[d a ,Z) = cfd ti , c t ? = d a Z» 

Then Ca are constants and 

CLL = C^ = 0. 

Furthermore 

[z,a] = - cz n 

where cz is either or 2. 

In addition, if Q is a null form, then 

(7.10) ZQ(dcj>, dtp) = Q(d<j), OZij) + Q(dZ(j), dij>) + Q(d<j>, ty) 
for some null form Q on the right hand-side. 

Proof. Since Z = Z a d a is a Killing or conformally Killing vector field we have 

(7.11) d a Zp + dpZ a = fm a p 

where Z a = m a pZ^ . In fact, for the vector fields above, / = unless Z = S in which case / = 2. In 
particular, 

L a L (3 d a Z f3 = 0. 
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If Ca is as defined above and c a p = Ccpm^p = d a Zp the above simply means that en = c— = 0. 
which proves the first part of the lemma. To verify (|7.1UI) we first consider the null form Q = Q a p We 
have 

ZQ af3 (d0, dip) = Q af s(dZ^ dip) + Qa^dcP, dZip) 

+ [z, daiwrf - d p <j>[z, d a ]ip + [z, ddWort - d a <p[z, d p ]ip 

= Q a( s{dZ4>, dip) + Qapidcj), dZip) - c^d^dpip - dp^ip) - (^(d^cpdaip - dafpd^ip) 

= Q a ?(dz<p, d^p) + Q aP {d<p, dZip) - c^Q^idcp, dip) - ^Q^{d<p, dip) 

The calculation for the null form Qo(d<p, dip) = m a ^ d a <pdpip proceeds as follows: 

ZQ (d<P, dip) = Q (dZ<P, dip) + Q (d<P, dZiP) + m a ?[Z, d^dpiP + m^d a <P[Z, dpty 
= Q (dZ<P, dip) + Q {d<P, dZiP) + m^c^d^dpiP + m^c^d^d^ 

= Q (dz& dip) + Q {d<j>, dZip) + fm af3 d a <pd p ip 

= Q (dZ(P, dip) + Q (d<p, dZiP) + fQ (d<p, dip), 

where / is a constant associated with a Killing (conf. Killing) vector field Z via a relation c a @ + c@ a = 
fm af3 . □ 

Lemma 7.5. If k 01 ^ is a symmetric tensor then 

(7.12) k a/3 [d a df3, Z] = kfd a dp, where kf = k a ^c^ + k^c", c/ = d a Z». 
In particular Iz^f = 2k a/3 and 

(7.13) \k z \ LL < 2\k\ LT . 
In general 

(7.14) [k aP d a d p ,Z I }= Y, k^dadpZ 1 ^ 

/i+/ 2 =/,|/ 2 |<|i-| 

where 

(7.15) k Ja ?= J2 c J K a ^ u Z K k^ = -Z J k a/3 - Y, Z K kf+ Y d^,Z K kT 

\K\<\J\ K+Z=J |if|<|J|-2 

for some constants ci?^ and d{?^ . Here the sum ()7.14j) means the sum over all possible order pre- 



serving partitions of the ordered multiindex I into two ordered multiindices I\ and Ii- 



Proof. First observe that since the vector fields Z are linear in t and x we have 

[dip, z \ = [9/3, Z]d a + [d a , Z]dp = c^d^d a + c^d^dp. 



which proves the first statement and the second follows since cfp = 0. 
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To prove (|7,14l) we first write 



K+J=I 



Then we observe that 



(7.16) Z J d a dp4> = I 2 ' 11 ' l ZH2 > [~> \ ZiXn -\ [Z Ll ",d 2 a(3 ]] ...] ] ] Z J *cP, 

Jl + J2 = J, Jl = (l>ll i-, l ln) 

where the sum is over all order preserving partitions of the ordered multiindex J = (lx, l^) into two 
ordered multiindices J\ = (in, t\ n ) and J2 = (^21, • ••> ^2fc)- It therefore follows that 

k JaP = - £ (Z K kr#) [Z«, [Z'\ [..., [z«-\[z\dl p ] ] ...] ] ] 

K+L=J,L={i 1 ,...,n) 



The desired representation follows after taking into account that 

{Z K k^)[Z,dl p ) = -{Z K kf)d a d p 



□ 



Corollary 7.6. Let U g = □ + H aP d a dp. Then with Z = Z + c z 
(7.17) U g Zcp - ZU g $ = -(ZH°P + Hf)d a l 

As a consequence, we have 

(7.W, |Ey* - zn a „ < (Ml + «) £ !«*'< 



|J|<1 



In general 



(7.19) □ 5 Z / 0-Z / D 5 = - ^ H^dadpZ 1 ^, 

h+h=i,\h\<\i\ 

where 

(7.20) ff jQ ^ = ^ clf^Z M H^ = -Z J H a < i - Z M Hf+ ^ d J ^Z M H^ 

\M\<\J\ M+Z=J |Af|<|J|-2 



We have 

(7.21) |D g ZV|<|i f 5 3 0[+ 1 + t 1 +|g| £ E l^ll^ 



|A'|<|/|, |J|+(|JS-|-1) + <|I| 

+ tt^..E.( E Jf J *i" + E i* J '»i-+ E 1^1)1^- 



\K\<\I\ \J\+{\K\-1)+<\I\ |J'| + (|A'|-l)+<|/t-l \J"\+{\K\-l) + <\I\-2 

where {\K\ - 1)+ = \K\ - 1 if |if | > 1 and {\K\ - 1) + = i/ |K| = 



30 



Proof. First observe that 

ZU g 4> = (Z + c z )U(P +{Z + c z )H a Pdlp4> 

= UZ<j> + H^dlpZcf) + (ZH^ffipt + (Hf + c z H^)d 2 ^ 
\p(j) + (Hz + c z ±i i" n << 



U g Z<t> + (ZH^)dl^ + [Hf + c z H^)d 2 



Recall now that the constant cz is different from only in the case of the scaling vector field S. 
Moreover, in that case 

Hf + c s H a P = 

The inequality (|7.18j) now follows from (|7,17|) , (|7,13|) and the estimate (|7.7|) . The general commutation 
formula (|7.19j) follows from the following calculation, similar to the one in Lemma [731 We have 

Z 1 ^ = Z J D0 + Z T H^dlp^ = UZ 1 ^ + Z J 'H°P Z k dlpcj) 

J+K=I 

If we now use ()7.16(l we get (j7.19|) as in the proof of Lemma 17.51 The inequality Q7.21JI now follows 
from (J7Hni), (j7T3l) and the estimate ((77)) . □ 

8 Basic energy identities 

We now establish basic energy identities for solutions of the equation 
(8.1) U g 4> = F 

We denote by T, t the hyper surfaces t =const, by C^(q) the forward light cones with a vertex at (q, 0) 
and truncated at times t\,t2. We also denote by K^(q) the interior of the light cone C 4 t 1 2 (g)and by B tjT 
the ball of radius r centered at (t, 0). 

Lemma 8.1. Let <f> be a solution of (|8.1j) . Then for any t\ < t2 and an arbitrary q < t2 
( - 5 00 |^| 2 + g^drfdrf) = [ ( - g 00 \M\ 2 + gVdiWjtf,) 
.2) L { d ^MM-\d t g^d a cpdpcp + Fd t( ^, 



and 



( - g 00 ^ 2 + gVdiWrf) + f t = [ ( - <?°°l^| 2 + fFdiWrf) 

+ 2 f (d a g a Pdp<t>d t <t> - \d t g aP d a 4>dp(p + Fdt 
(8.3) +2 I {2(g<*P - m^)L a dp^ + {g af3 - m^)d, 
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Proof. We multiply the equation (|8,1|) by dt4> and integrate over the space-time slab between the hyper 
surfaces S tl and S i2 . We have 

- r / S^Sfatdtt = f 2 f {g a(3 d^d t d a ^ + d a g^d^d t 4>) 

g°%<l>d t <l> + f g^dpW 
\f ( - 3°°l^| 2 + g'Uho^o) -\\ ( - 5 00 |a^| 2 + g^di^cf) 



+ 



( 2 ! {da^dpMf - \d t g^d a 
and the desired identity ()8.2|) follows. Similarly, integrating over the region K^(q) we obtain 
g^d^ + g^d^djcf) - [ t _ {2g^L a d p <pd t <t ) + g^d^dp ( t>) 

{-g 00 \d^\ 2 +g^d t ^) 



+ 2 



< 2 (<7) 



Subtracting the Minkowski part from the metric g in the CS(q) integral leads to the identity l|8.3|) . □ 
Corollary 8.2. Let 4> be a solution of the equation (|8,1|) with a metric g satisfying the condition that 

(8.4) \ H \^J> H a/3 = g a/3 - m a/3 . 

Then for any < 7 < 1 

r*2 r ~ l/U|2 



(8.5) + 8 [ 2 [ \d a g a(S d^d t 4> - \d t g a ^d a (pd^ + Fd t <f>\ 



Proof. First we note that Q8.4|) implies that 

(8.6) ^(|d t( />| 2 + |V,0| 2 ) < -<? 00 |^| 2 +g ij d i <l>d j <j> < -({d^ 2 + |V^| 2 ) 
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The inequalities (|8,3j) and (|8. 2 (l imply that 



(8.7) 



/ \8<P\ 2 <[ (-g^dt^ + g^d^) 

Jcll{q) JZu 



+ 2 I t d a g^d P 4>d t 4> - \d tg ^d c 
'< 2 W 2 



+ Fd t 



< 



+ 2 



+ 2 / 2tf* - m ap )L a d^d t( l> + {g a<i - m^)d, 

gW^+g^d^) 

d a g^d^dt^ ~ \d t g a/3 d a 0d^ + Fd t <t> 
+ 2 / 1 2(g a ? - m a P)L a dp<t>dt<t> + {g aP - m^)d, 

JC&lq) 



t2 



t T(?) 

We multiply the above inequality by an integrable factor 7(1 + \q\)~ l ~ 21 and integrate with respect to 
q in the interval (—00, to obtain: 



*2 r 
h Js t ( 



(8.10) 



+ 2 T [ \d a g° 
Jti JE r 

+ 2 Z La + k\) 1+2 ^ 



\d t g a<i d a c\>d^ + F9^| 



7 



_ m^d^d^ + 2{g^ - m^d^dtcP 



LP _ ^LP\ 



where we also used (j8.6[) . On the other hand using (|8.6I) and (|8.2j) yields 

(8.11) / {m 2 + |v^| 2 ) <\l (\d t( p\ 2 + \v x 4>\ 2 ) 



.12) 



and the corollary follows. 



+ 



ai3 



\d t g aP d c 



+ F8 t 



□ 



9 Poincare and Klainerman-Sobolev inequalities 

We now state the following useful version of the Poincare inequality. 

Lemma 9.1. Let f be a smooth function. Then for any 7 > —1/2,7 1/2 and any positive t 

\f(x)\ 2 d X ^ f . f . 2jo , ^ f \d r f(x)\ 2 dx 



(9.1) 



+ |t-r|) 2 + 2 7 
R 3 S(t+i) 



<c J \f\ 2 dS + C J ± 

S( t +i) R 3 



+ \t- rj) 2 ^ 



provided that the left hand side is bounded. Here Sr t+i ) is the sphere of radius t + 1 



ana r = x . 
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Proof. Using polar coordinates x = rto we write 

rt+l 



rt+i 

\f(r,oj)\ 2 -\f(t + l,Lo)\ 2 = -2 / arf(p,u)-f(j>,u)dp 

J r 



Hence 

rt+l 



ft+1 

\f(r,Lo)\ 2 r 2 <\f(t+l,u;)\ 2 (t + l) 2 + 2 / \d r f(p,u)\ \f(p,u)\ p 2 dp, if r<t+l. 

J r 

Therefore multiplying by (1 + \t — r|) -2-27 and integrating with respect to r from to t + 1: 



f t+1 \f(r,u)\ 2 r 2 dr < + 1, u;)| 2 (t + l) 2 dr f t+1 f t+1 \d r f(p,u)\\f(p,u)\ 

J (l + \ t -r\) 2 + 2 ^ J (l + | t - r |)2+27 + J Q I 



(l + \t-r\) 2 + 2 ~< 

<i/(t+i, W )i (t+i) +y o y o (1 + |f _ r|)2+27 dp 

< + 1, ,)| 2 (t + i) 2 + f 1 p 2 

<im i i i) 2 i ( r 1 \Mk^)ll± Y /2 ( f t+1 \f(p^)\ 2 p 2 dp y/ 2 



where we first changed the order of integration and then used Cauchy-Schwarz inequality. It therefore 
follows that 

\f(r,u)\ 2 r 2 dr \d r f(p,u J )\ 2 p 2 dp 

and if we also integrate over the angular variables we get 

f \f(x)\ 2 dx < /■ 2 /■ |d r /Qr)| 2 rf:r 

7 (l + |t-r|) 2 + 2 7 ~ 7 l/|a " + y (l + |t_ r |)27 

|x|<(t+l) S (t+1) \x\<(t+l) 



On the other hand, if we instead integrate from t + 1 to 2(i + 1) we similarly obtain 
/■W) )/(r^)l 2 r 2 rir < Z^ 1 ) + 1, a;)| 2 (t + l) 2 ^(t+i) yr |g r/(j0 ^)[| /(;0 , 

y t+ l (l + |t-r|) 2 + 2 ^y t+1 " (l + |t_ r |)2 + 2 7 - + J w J w (1 + | t _ r |)2 + 2 7 

2(t+l) f2(t+l) 



sl/(t+1 ,^ (t+1 , + £-' i^y ^ 



and as before it follows that 

\f(x)\ 2 dx . f |J|2jc , /■ |<9 r /(x)| 2 dx 



(l + \t-r\)^ 

(t+l)<|x|<2(t+l) S (t+ i) (t+l)<|x|<2(t+l) 
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Finally, in the region r > 2(t + 1) the estimate (|9.1|) would follow from the Hardy type inequality: 

<M> / / «^ + ( t+ l)-^ / 1/l'dS. 

Nl>(t+i) M>(m) 

that hold provided the left hand side is bounded. One can for the proof assume that / ha compact 
support since we can choose a sequence of compactly supported functions converging to a given function 
/ in the norm defined by the right hand side as long as the norm in the left of / is bounded. (|9.2[) for 
compactly supported smooth functions can be easily seen from integrating the identity 

/ 2 f 2 \ 2 r 2 r 2 
°r {^k) = ■ W + (1 " 2 7)-^/ 2 > 7 + "1/2 

from r = t + ltor = oo and using Cauchy-Schwarz as above. □ 

We now state the global Sobolev inequality, which is due to S. Klainerman |Klj . 

Proposition 9.2. The following inequality holds for an arbitrary smooth function <j>. 

\</>(t, x )\(l + t+\t- r|)(l + \t- r\) 1 ' 2 < C ^ WZ'Ht, OHtf- 

|/|<3 

10 Decay estimates for the wave equation on a curved space time 

In this section we will derive some basic estimates for the scalar wave equation on a curved background. 
The results will require some weak assumptions on the metric g, which will be easily verified in the 
case of a metric satisfying the reduced Einstein equations. 
We consider the reduced scalar wave equation: 

(10.1) UgCj) = F. 

The following result is a generalization of the lemma in |Llj to the variable coefficient case: 

Lemma 10.1. Suppose that 4> satisfies the reduced scalar wave equation (jlO.lj) on a curved background 
with a metric g. Suppose that H a @ = g a @ — m a ^ satisfies 

(10.2) \H\ < ~, and \H\ LT < \: ' <7 ' + ' 



4' 1 1 41 + i + \x\ 

when t/2 < \x\ < 2t and 

f°° dt 1 

(10.3) || Hit, .)\\ Loo{Dt) — t < -, where D t = {x € R 3 ; t/2 < \x\ < 2t}. 

Then for any t > and x£l 3 , 

(10.4) (l + t+\x\)\d<t>(t,x)\<C sup Y \\Z^(t, -)\\ L oo 

0<r<t 



+ C f ((1 + t)\\F(t, OIU-^) + ^ (1 + r^H^r, Olk-CD.)) dr 

J ° |J|<2 
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Proof. Since by Lemma 17.21 

(10.5) (l + \t-r\)\d(j)\ + (l + t + r)\d(j)\<Cj2 \ zI( t>l r = \x\, 

\i\=i 

the inequality (|1U.4|) holds when r < t/2 + 1/2 or r > 2t — 1. Furthermore, since 

(10.6) (l + r)\d q 4\<C\d q {r<l>)\ + C\<l>\, r>l 
it follows that 

(10.7) (1 +t + r)\d(f>\ < C \Z T <t>\ + C\d q {r<f>)\ 

|/|<i 

Hence it suffices to prove that \d q {rcj))\ is bounded by the right hand side of ()10.4j) when t/2 + 1/2 < 
r < 2t — 1. By Lemma 17.31 



lOi 



5 l"l ' |/|<2 



and using the decay assumptions (|10.2|) and (|10.7|) we get 
(10.9) 

\(4d a - §^d q )d q (r<f>)\ < JfL|9 ? (r0)| + ^ + C(i + 1)|F|, when t/2 + 1/2 < r < 2t - 1 

5 |/|<2 

Along an integral curve (t,x(t)) of the vector field d s + H—(2g L —)~ 1 d q , contained in the region 
t/2 + 1/2 < |x| < 2t — 1, we have the following equation for ^ = d q (r<p): 

(10.10) <%>| + / 

where h = C\H\/(l + 1) and / = Ct|F| + CJ2\i\< 2 l^l/C 1 + 0- Hence multiplying (tTinUl) with the 
integrating factor e~ H , where H = f h(s) ds we get 



(io.il) f t ^ e ~ H ) - fe 



-H 



If we integrate backwards along an integral curve from any point (t,x) in the set i/2+1/2 < |x| < 2t— 1 
until the first time the curve intersects the boundary of the set at (r,y), \y\ = t/2 + 1/2 or | y \ = 2 r — 1 , 
we obtain 

m,x)\ < exp [|fc((7, -)|U- dff) y)l + y exp ||/*(a, dtr) ||/(r', dr', 

where the L°° norms are taken only over the set t/1 + 1/2 < |x| < 2t — 1. (Note that any integral curve 
has to intersect either of the two boundaries r = t/2 + 1/2 or r = 2t — 1 since the slope of the curve x(t) 
has to be close to 1 when Hll is small.) The lemma now follows from taking the supremum over x in 
the set t/2 + 1/2 < |x| < 2t — 1, using that on the cones \y\ = t/2 + 1/2 or \y\ = 2t — 1 we have that 
M < Cr\d q 4>\ +C\4>\ < CJ2\i\<i \ zI <f>l b y and using that by (fTIHl) Jj ||/i(cj, -)||l- At < |. □ 
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For second order derivatives we have an estimate which gives a slightly worse decay: 



Lemma 10.2. Let (ft be a solution of the reduced scalar wave equation on a curved background with a 
metric g. Assume that H a ^ = g a ^ — m a/3 satisfies 



;io.i2) 



E \ zlH \ 

\i\<i 



< 



and 



UI<1 



+ t+ X 



when t/2 < \x\ < It for some i < 1. Then, for t > 0, x £ M 3 , we have 



(10.13) (l + i + |x|)V |SZ 7 (/>(t,x)| < C sup (- 

— ' f1<-r</ ^ 1 

I/I<1 

"* / 1 + t\ GS 



l + t\Ce 



+ T 



Ell^>0IU° 

|/|<2 



c 



1 + T 



E (1 + r)||^F(r, 0IU~ (JJt) + + O^H^^ OIU-p.)) <*r, 

|/|<1 |/|<3 



where D t = {x G R 3 ; t/2 < |x| < 2t}. 

Proof. First when r < t/2 or r > t/2 the lemma trivially follows from (|1U.5|) with (ft replaced by Z<f> so 
it only remains to prove the lemma when t/2 < r < 2t. We have 

(10.14) UgZcj) = F z = ZF + {U g Z<t> - ZUg<j>) , 

where by (|7.18|) the additional commutator term can be estimated by 



;io.i5) \u g z<t>-zn g <i>\< 



ZH\ + \H\ , \ZH\ LL + \H\ LT 



l + t + \q\ 



+ 



y \ dz ^\ < — £ - — y \dz\ 

|/|<l |/|<l 



where we used the decay assumption (jl().12|) . Furthermore with the help of (|l().7jl . applied to Z (j) in 
place of 4>, we obtain 



;i0.16) 



\n g z<j>-zn g 4>\<- 



+ t + 



E 1^(^)1 + E i^O' 

|7|<1 |J|<2 



Hence by (fTTH?)) applied to (|1U.14|) in place of (fT0~Tj) we get 



(10.17) 



1(45, 



H_LL 

25 



LL 



d q )d q (rZcf>)\< E V 



J <3 



+ t 



_£|0 fl (rZ'fl|+i(|ZI?'| + |*'|) 
|i|<i 



when t/2 + 1/2 < r < 2t - 1. Therefore 



(io.i8) | (40. - e i»«(rzV)i| < Yh E + c E rrr + ct(|ZjF| + 

|/|<1 |J|<1 |/|<3 

The desired result follows multiplying ()10.17|) by the factor (1 + i)~ c ^ and integrating as in the proof 
of the previous lemma. Along an integral curve we have the equation 



10.19) 



d 



^(v(i+r C£ ~) <(i+ty u£ f, 



-Ce. 
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where 

(10.20) V=E / = C(l + t)(|^| + |if|) + CX;^ 

l-f|<l |J|<3 

The lemma now follows as in the proof of Lemma 110.11 □ 
We observe that similar estimates hold for a system 

(10.21) Ogfa = 

In particular, in our case, certain components of Fu U expressed in the null-frame will decay better than 
others and for these components we will also get better estimates for <j)^ v . Since the vector fields L and 
L commute with contractions of any of the vector fields {L,L, Si, S2} proofs of the preceding lemmas 
imply the following result: 

Corollary 10.3. Let <fiuv be a solution of reduced wave equation system (jlU.21f) on a curved background 
with a metric g. Assume that H a @ = g a @ — m a ^ satisfies 

(10.22) £i* z #i<§, and E \ zlR \^ + i g br ± \ rrfeid - 

i/i<i i/i<i 1 1 

when t/2 < \x\ < It, for some e < 1 and 

r°° dt e 

(10-23) ^ ||^.)|| LOO(DO __<_. 

where D t = {x G R 3 ; t/2 < \x\ < 2t}. Then for any U,V G {L,L, Si, S 2 } and any t > 0, x G K 3 : 
(10.24) 

(l+t+\x\) \d(f>(t,x)\ uv < C sup V ||Z 7 0(t,-)IU- 

0<r<* |J| < 1 

+ c I [{i + t)\\\f\ uv (t,.)\\ lx{Dt) + ^ (i + t)- 1 !^ 7 ^,-)!!^^))^, 

|/|<2 

(10.25) 

(l + t + |x|)^ < C sup (i±l) C£ ^||Z^(r,.)||L- 

l-f|<l l-f|<3 

Proof. By Lemma 17.31 for each component we have the estimate 



(10 , 6 ) I («, - |!H - ^|±fi) 8 ,(^) + gs| < (1 + + wy> £ 1*' 



/«/ I 

I\<2 
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and since d s and d q commute with contractions with the frame vectors L,Lwe get 

(io,7) |(4a, - ^ Bq - ^±^) dMuv) + ^\< {l + r_m^ + m y> E 



|J|<2 



As before it also follows that 

(10.28) (1 + t + \r)\dcl>\ uv < \ Z ^\ + \9 q (r4>)\uv 

|J|<1 

The lemma now follows as before. □ 

11 Energy estimates for the wave equation on a curved space time 

In this section we derive the energy estimate for a solution cf) of the inhomogeneous wave equation 

(11.1) Ug4> = F 

under the following assumptions on the metric g a/3 = m a ^ + H a/3 : 

(1 + \q\y l \H\ LL + \8H\ll + \8H\ < Ce(l + t)~\ 

(11.2) (1 + {qiy 1 \H\ + \8H\ < Ce(l + t)-*(l + \q\)~^ 

Proposition 11.1. Let <p be a solution of the wave equation (|11.1() with the metric g verifying the 
assumptions 1)11.2(1 . Then for any < 7 < 1/2, there is an eq such that for e < Eq 

Remark 11.2. Observe that by the Gronwall inequality the energy estimate of the above proposition 
implies t £ growth of the energy. For similar estimates, proved under different assumptions, see also 

Proof. The proof of the proposition relies on the energy estimate obtained in Corollary 18.21 

'■t r „, 



+ 8 / f \d a9 ^d^d t <p - htg^d^d^ + Fd t 

Jo Jt, t 1 
t 



2 I I M Z MU+2-y - m a P)d a <t>dp<t> + 2(g^ - m^)fy<W| 
Jo Jt, t K 1 + \Q\) ' 

We start by decomposing the terms on the right hand side with respect to the null frame. 

{d^d^d^l < (\H\ \8H\ + \(8H)ll\ + \BH\) m 2 + \8H\ \B<t>\ |d^| 

Similarly, 

Idtg^d^d^l < (\g - m\ \8g\ + \{dg) LL \ + \8g\) |50| 2 + \dg\ \B<t>\ \d<j>\ 
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Therefore, using the assumptions (|11.2|) on the metric g, we obtain that 

(n.4) \d a g^d P W - \d t g ap d a <pd^\ < j^l^l 2 + (1 + [ ^ |)1+27 l^l 2 

Decomposing the remaining terms we infer that 

\(g a/3 - m^daWptl < \H LL \ \dcp\ 2 + \H\\5<I>\ |d<£| 

Similarly, 

_ m ^) La d^d t 4>\ < \h ll \ m 2 + \h\ m m 

Once again, using the assumptions (|11.2|) . we have 

(11.5) \2{jf# - m^)L Q d p( f>d t <p + {g^ - m^)d a <j>dp<j>\ < e^f\dch\ 2 + - £ \5<p\ 2 

Thus 

k, JO k r (1 + l«l) 1+27 " J* k kM+t (1 + l«l) 1+21 ^ JO k, 

and the desired estimate follows if we take e so small that Ce < 7/2. □ 

12 Estimates from the wave coordinate condition 

In previous sections we have shown that one only needs to control certain components of the metric in 
order to establish decay estimates for solutions of the reduced wave equation. In this section we will 
see that the wave coordinate condition allows one to estimate precisely those components in terms of 
tangential derivatives or higher order terms with better decay better. Recall that the wave coordinate 
condition can be written in the form 

(12.1) ^(<r^det^|)=0 
We have the following decomposition 

g^y/\detg\ = [mT + H» v ) (l - hssH + 0{H 2 )) 

where H a/3 = g a/3 — m a/3 , h a g = g a p — m a R. Recall also that g a @ is the inverse of g a p and H a @ = 
—m fJ-a m 1/ P + 0(h 2 ) . Therefore we obtain the following expression for the wave coordinate condition: 

(12.2) d^H^--m^ti H + 0» v (H 2 )} =0 

Using that we can express the divergence in terms of the null frame 

(12.3) = L^d q F» - L^d s F» + A^aF" 
we obtain: 
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Lemma 12.1. Assume that \H\ < 1/4. Then 

(12.4) \8H\ LT <\dH\ + \H\\dH\, \dtrH\ < \&H\ + \H\ \dH\ 
Proof. It follows from (TEH?]) and (TEH]) that 

(12.5) \L^d(H^ - -m^ tr H) \ < \dH\ + \H\ \dH\ 

Contracting with T u and using that rriTL = gives the first inequality and contracting with and 
using that ttill = — 2 gives the second since 

(12.6) H LL + tr H = \xH 

□ 

We now compute the commutators of the wave coordinate condition with the vector fields Z. 

Lemma 12.2. Let Z be one of the Minkowski Killing or conformally Killing vector fields and let tensor 
H satisfy the wave coordinate condition. Then the estimate 

\dH J \ LT < Yj \dZ J H\+ \Z h H\ ■ ■ ■ \Z h H\ \dZ h H\ 

\J\<\I\ h+...+I k =I,k>2 

holds true for the expression 

(12.7) H^ = Z J H^+Y^ (/j^Z J H^, where H^ = H^-hn^trH 

\J\<V\ 

with some constant tensors cf 111 such that Cjf = if \ J\ = \I\ — 1. 

Proof. The wave coordinate condition ([12.1)1 can be written in the form 

d„(& v ) = 0, where = {m^ + H^)^/\detg\. 

Let Z be one of the Minkowski Killing or conformally Killing vector fields. Then for any vector field 
F we have that 

Z I d a F a = d a (z I F a + cj«Z J F^=d a ( £ cj«Z J F^, 

\J\<\i\ \J<\i\ 

where cj^ are constants such that 

c j 7 « = ( r 7 , if \J\ = \I\ and cj L - = 0, if \J\ = \I\-1 

The last identity is a consequence of the relation between cj a 7 and the commutator constants c Q/ g = 
[d a , Z]p for which we have established that cll = 0. It therefore follows that 

M E 4"Z J Gv,)=0. 

\J\<\I\ 
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Decomposing relative to the null frame (L, L, Si, S2) we obtain 

|J|<|J| |J|<|/| ' |J|<|/| 

We now contract the above identity with one of the tangential vector fields T v , T € {L, S\,S2} to 
obtain 

V'T^dgZ'G^+Y, Cj^V 8 q Z J G lu \< J2 Wg\ 
\A<V\ \A<V\ 

We examine the expression 

L 1 T v Z J d q G^ v = V T v d q Z J ({m lv + H^)J\&Ag\ \ = £ V T v 3 q ((Z h H lu )Z j2 y/\deti\\ 

\ ' .h+j 2 =j \ ' 

since ttilt = T u m-y V = 0. The desired estimate now follows from the identity \f \ det g\ = 1 + f(H), 
which holds with a smooth function f(H) such that f(H) = —trH/2 + 0(H 2 

We now summarize the results of this section in the following 

Lemma 12.3. For a tensor H obeying the wave coordinate condition 

(12.8) \dH\ LT < \8H\ + \H\ \dH\, 

and 

(12.9) 



□ 



\8ZH\ LL < \dH\ LT + E \9Z Z H\ + \ Z * H \ \ dzJR \ 

\i\<\ \i\+\J\<y 



In general, 



(12.10) \dZ l H\ Lr < \dZ J H\ + £ \dZ J H\ + £ \Z Im H\ 

\J\<\I\ \J\<\I\~1 |/i|+...+|/ m |<|/|,m>2 

and 



\Z h H\\dZ h H\, 



(12.11) |dZ 7 iI| LL < £ \dZ J H\+ £ \9Z J H\ LT + Yl \ dzJR \ 

\J\<\I\ |J|<|/|-1 l^l<l-f|-2 

+ \Z Im H\- ■ ■\Z h H\\dZ Il H\. 

|/i|+...+|/ m |<|/|,m>2 

The same estimates also hold for H replaced by h. 

Proof. This follows directly by the previous lemma with the help of the identities tolt = and 
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13 Estimates for the inhomogeneous terms 

In this section we will show that the inhomogeneous terms of the reduced Einstein equations can be 
estimated in terms of tangential derivatives, for which we have better decay estimates, or tangential 
components which in turn can be expressed, using the wave coordinate condition, in terms of tangential 
derivatives and lower order terms. Recall that according to Lemma 13.21 the symmetric two tensor 
hfiu = 9[iu — fn^v verifies the reduced Einstein equations of the form: 

□s^V = Ffiv{h)(dh,dh), 

(13.1) F^{h){dh,dh) =p(dph,d v h) + Q^(ah,dh) + G tu ,(h)(dh,dh), 

(13.2) P(d^k, 8 uP ) = ^tvk d v tvp - ^k^d^, 

Here Q\iv are linear combinations of the null- forms and G/ M/ (h){dh, dh) is a quadratic form in dh with 
coefficients that are smooth functions of h and vanishing at h = 0. 

Lemma 13.1. The quadratic form P satisfies the following pointwise estimate: 

(13.3) \P(dp,dk)\ ru < \dp\\dk\ + \dp\\dk\, 

(13.4) \P(dp,dk)\ < \dp\ TU \dk\ru + \dp\LL\dk\ + \dp\\dk\ LL 

Proof. The first part of the statement follows trivially from (|13.2|) . To prove ()13.4j) we use (|5,9|) applied 
to R^d^p in place of p and S v d v k in place of k, for any vector fields T and S, to obtain 

(13.5) \T»S V P{d^p,d v k)\ < \T^d^p\ TU \S u d v k\ru + \T^p\ LL \S u d v k\ + |T^| \S v d v k\ LL 
which proves the lemma. □ 



Using the additional estimates on the hn component, derived in Lemma [1231 under the assumption 
that the wave coordinate condition holds, we obtain the following: 

Corollary 13.2. Under the additional assumption that h satisfies the wave coordinate condition (|3.4|) . 
the quadratic form P obeys the estimate 

(13.6) \P(dh,dh)\ T u<\dh\\dh\, 

(13.7) \P(dh, dh)\ < \dh\\ u + \8h\ \dh\ + \h\ \dh\ 2 

Moreover, 

|Z 7 P(a/i,a/i)| < {\dZ J h\ TU \dZ K h\ ru + \dZ J h\\dZ K h\) + \dZ J h\ LT \dZ K h\ 

\J\+\K\<\I\ I-/| + |A'|<|/|-1 

+ Yj \dZ J h\\dZ K h\+ Y \Z Jm h\- ■ -\Z- h h\\dZ j2 h\\dZ Jl h\ 

|J|+|Jr|<|Z|-2 \Ji\+... + \J m \<\I\,m>3 

Proof The inequality (|13.6|) follows directly from (|13.3|) . To prove ()13.7(l we use ()13.4|) and that by 
the wave coordinate condition \3H\ll < \dh\ + \h\ \dh\. 
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We now note that Z 1 P(d^h,d v h) is a sum of terms of the form P{d a Z J h : dpZ K h) for some a, j3 
and | J| + \K\ < I: 

iZ^idh^dh)] <C Yj \P{dZ J h,dZ K h)\ 

\J\ + \K\<\I\ 

It follows from Lemma 15.41 and Lemma 112.31 that 

(13.8) \P(dZ J h,dZ K h)\< \dZ J h\ ru \dZ K h\ru + \dZ J h\ LL \dZ K h\ 

\J\ + \K\<\I\ \J\+\K\<\I\ 

< Yl \ Bz ' h \ \ zK]l \ + \ dz ' 1h \ru \dZ K h\ ru 

\J\ + \K\<\I\ 

+ zZ ( 12 \ dzJ ' h \cr + Y \ dzJ " h \ + Yl \Z Jm h\ ■ ■ -\Z- h h\\dZ- h h\}\dZ K h\ 

|J1+|-K"l<|i"| |J'|<|J|-1 \J"\<\J\-2 \Ji\+... + \Jm\<\J\,m>2 

which proves the lemma. □ 

Proposition 13.3. Let F^ u = F fJiU (h)(dh, dh) be as in Lemma \3. "A and assume that the wave coordinate 
condition holds. Then 

(13.9) \F\ru < \dh\ \dh\ + \h\ \dh\ 2 
and 

(13.10) \f\ < \dh\ 2 TU + \8h\ \dh\ + \h\ \dh\ 2 



(13.11) \ZF\ < (\dh\ru + \dh\ + \h\ \dh\) (\dZh\ + \dh\) + (\dZh\ + \Zh\ \dh\) \dh\ 

(13.12) < Y {\dZ J h\ ru \dZ K h\ru + \dZ J h\\dZ K h\) + Y \dZ J h\ LT \dZ K h\ 

\j\+\k\<\i\ |j]+|i<r]<|/|-i 

+ Y \dZ J h\\dZ K h\+ Y \Z Jm h\- ■ ■ \Z j3 h\ \dZ j2 h\\dZ Jl h\ 

\J\ + \K\<\I\-2 |Jx|+-+|Jm|<|/|,m>3 

Proof. First 

|^ 7 G^(/i)(a/i,a/i)| < c Y \ zhh \ • ■ • \ zhh \ \dZ h h\ \dZ h h\. 

|/i|+...+|4|<|/|,fc>3 

Since ZQ(du,dv) = Q{du,dZv)+Q(dZu,dv)+a^Q ij (du,dv), and \Q^{dh,dk)\ < \dh\ \dk\ + \dk\ \dh\ 
it follows that 

|Z 7 Q^(a/i,a/i)| < C Y \Q^(dZ J h,dZ K h)\ < C Y \dZ J h\\dZ K h\ 
\J\+\k\<\i\ \A+\k\<\i\ 

□ 
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14 The decay estimates for Einstein's equations 

In this section we will establish the improved decay estimates for Einstein's equations. Our strategy 
is to use the weak decay estimates, obtained from the assumed energy bounds, to prove sharper decay 
estimates and then to recover the energy bounds in the next section. 

Theorem 14.1. Suppose that for some < 7 < 1/2 



(14.1) \dZ J h\ < Ce(l + t+ |g|)- 1/2 - 7 (l + |q|)" 1/2 " 7 , |/| < N/2 + 4, 

(14.2) \Z'h\ ^Ceil + ty 1 , q = l, \I\ < N/2 + 4 

hold for < t < T. Then for < t <T we have 

(14.3) \Z l h\ < Ce(l + t + |q|)~ 1/2 - 7 (1 + k|) 1/2 " 7 , |J|<iV/2 + 4, 

(14.4) \dZ J h\ <Ce(l + t+\q\)- 3/2 ^{l + \q\) 1/2 ^ f , \I\ < N/2 + 3. 



Assume also that h satisfies the wave coordinate condition. Then for < t < T we have 

(14.5) \dh\ LT + \dZh\ LL < Ceil + t)- 1 - 2 ^, and \h\ LT + \Zh\ LL < Ce(l + 1)" 1 ^ + |g|). 

Furthermore if in addition h satisfies Einstein's equations then for e sufficiently small and < t < T 
we also have 



(14.6) \dh\ru <Ce(l + t)-\ \ h \ TU < Ce(l + t)~\l + \q\), 

(14.7) < C^(l + 1) -1 ln(2 + 1). 

In general, there are constants Mk, Ck and Ek > such that if e < St then for \I\ = k < N/2 + 2 

(14.8) IdZ 1 ^ < C k s(l + t)- 1+Mfc£ , and \Z J h\ < C k e{\ + t)- 1+M " £ {l + \q\) 



Remark 14.2. We remind the reader that, as stated in the Remark I2.4[ our estimates make no dis- 
tinction between the tensors h and H = —h + 0(h 2 ). In particular, one can directly verify that the 
conclusions of the theorem also hold for the tensor H. 

First we note that all the estimates (|14.3j) - (|14.8j) trivially follow from the assumptions (|14.1|) - (|14.2[) 
away from the light cone, thus the theorem is only useful in the region t/2 < \x\ < 2t. The estimate 
([14.3ft follows from integrating (|14.1|) from q = 1, where ()14.2j) hold. Similarly the second parts of 
(1T43|) . (fTH)|l and (jT4"%|) follow from integrating the first and using (fTT^ . It follows from (|TO|) 
and Lemma 17. 21 that we have the better estimate (|14.4[) for the derivatives tangential to the outgoing 
Minkowski cones. The inequalities (|14,5[) - (jl4.8j) for tangential derivatives certainly follow from ()14.4j) . 
so it only remains to prove these estimates for a derivative transversal to the light cone. 
The missing improved estimates for a (dt — d r ) derivative transversal to the light cones will be obtained, 
in the case of (jl4.5j) . from the wave coordinate condition, see section WH and for H14.6|) - ([14.8j) . from 
integrating the reduced Einstein wave equations, see section llOl The estimates from the wave coordinate 
condition are easily obtained. In fact the first estimate in (|14.5|) follows directly from Lemma 112.11 
using the estimates ([14. lj) . (|14.3[) and ([14. 4j) and the second estimate in (|14.5[) follows integrating the 
first from q = 1 where ([14.2J) holds. However, the wave coordinate condition does not give estimates 
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for a transversal derivative of all components of the metric and the remaining components have to 
be controlled by integrating the wave equation expressed in polar coordinates. The estimates for the 
transversal derivative obtained from the wave coordinate condition rely on a decomposition of the 
metric with respect to the null frame. On the other hand, the estimates obtained from integrating the 
wave equation are based on a decomposition of the wave operator in terms of tangential derivatives 
and a transversal derivative. 



14.1 Proof of (fTOjl and rfTOjl 



For a fixed angular variable uj we integrate in the radial direction and use 1)14.1(1 and 1)14.2(1 

H+l 



rt+i 

14.9) \Z I h(t,r,u))\<\Z I h{t,t + l,u)\+ I |a r Z J /i(t, p,u)\dp 

J r 



< f t+1 Cedp < _Ce_ + Ce(l + \t - r]) 1 / 2 ^ 



1 + t J r (1 + t + \t - /9|) 1 /2+T(l + \t - p|)l/2+7 ~ 1 + t (1 + t + r)V2+7 

The estimate ()14.3() now follows. By Lemma 17^1 and ()14.3() 

which proves 1)14.4)1 . 
14.2 Proof of (11431) . 

We now show that the wave coordinate condition allows one to control certain components by lower 
order terms and terms with fast decay. 

Lemma 14.3. Suppose that the estimates (|14.1j) - (j!4.4|) hold and that h satisfies the wave coordinate 
condition. Then 

(14.10) J2\ dZlh \ LL+ E \ dZ ' h \cr< \dZ K h\+e{l + t + \ q \)- 1 -'^ 

\I\<k |J]<jfc-l \K\<k-2 

Here the sum over k — 2 is absent if k < 1 and the sum over k — 1 is absent ifk = 0. Furthermore 
(14.11) 

z J—(Y,\Z I h\ LL + J2 \Z J h\cr+ Y \Z K h\)(t,x)< sup \dZ K h(t,y)\+-^- 

L + lQl \I\<k \J\<k-l \K\<k-2 t/2<\y\<3t/2 imk _ 2 l+« 

Proof. We first prove (|14.1Uj) . Using the estimates of Lemma 112.21 derived from the wave coordinate 
condition followed by (|14.1[) - 1)14.4)) we obtain 

(14.12) 

Y |<9Z 7 F| LL + Y \dZ J H\ LT < Y \dZ r h\+ \dZ K h\+ Yl \Z Im h\-\Z^h\\dZ h h\ 

\I\<k \J\<k-l \I\<k \K\<k-2 |7i|+...+|/ TO |<fc,m>2 

< Y \dZ K h\ + e(l + t+\q\)- l ~ 2 ~< + e{l + t+\q\)- 1 - 2 ^ 

\K\<k-2 
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The proof of estimate (|14.11|) for |g| > t/2 follows directly from (|14.3j) . Thus we may assume that 
\q\ < t/2. We now use the inequality 

(14.13) \H(t,ru)\<\H(t,(t + i)w)\ + (l + \q\) sup \d p H{t, {t + p)u)\, 

|p|<kl+i 

and the boundary condition (|14.2f) to conclude that 

(14.14) \Z I H\LL + \Z J H\L r + \Z K H\ ii ^ z i H \ LL + ^ z J H \ LT + ^ z K H ^ y) + ,_ c 



1 + l<?| t/2<\y\<2t V / L + t 

The desired result now follows from (|14.10|) . □ 

The first part of 1)14. 5j) now follows directly from the Lemma with k = 0, 1 and the second part 
follows from integrating the first and using the boundary assumption (|14.2|) as in the proof of (|14.3|) . 

14.3 Proof of (fT4l)j) - (fT4~7jl . 

We will appeal to the L°° estimates of section^] for the reduced wave equation 

where F^ u is as in Lemma 13.21 We will now prove (j!4.6f) and (j!4.7f) assuming (|14.1|) - (|14.5|) . 

Lemma 14.4. Suppose that the assumptions of Proposition 1\ hold and let F^ v = F^ u (h)(dh,dh) 
be as in LemmaUTM Then 



(14.15) \F\ru < Cet~ l - 2 ~<\dh\ 
and 

(14.16) |F| < Cer l - 2l \dh\ + C\dh\ 2 ru 

Proof. This follows from Lemma 113.31 using (j!4. 1(1 - 1)14. 51) . □ 

Using the first part of Corollary 110.31 (|1U.24|) . and (|14. I|) - (jl4.5|) and the previous lemma we get 
Lemma 14.5. With a constant depending on 7 > we have 

(14.17) (l + t)\\\dh\Tu(t,-)\\ L ™ <Ce + Ce f\l + T)-^\\dh(r, dr, 

Jo 

and 

(14.18) (1 + t)\\dh(t, <Ce + C J* (e(l + T)-^\\dh(r, + (1 + r)\\\dh\ TU (T , -)\\ 2 L J) dr. 

The estimates ()14.6|) and (|14.7|) now follow from the above lemma and the following technical result 
applied to noo(i) = (1 + *)|| \dh\<ru(t, Oik- and ™oi(*) = (l + t)||0fr(i,-)lk<»: 



47 



Lemma 14.6. Suppose that n 00 > and n i > satisfy 

(14.19) n 00 (t) < Ce (J o (l + s) _1 ~ 7 ^oi(s)ds + l) 

(14.20) noi(i) < Ce (J Q (l + s) _1_7 n i(s)ds + l) (1 + s)- 1 n 00 (s) 2 ds 

for some positive constants such that < 16(C 2 + C)e < 7 < 1. TTien 

(14.21) n 00 (t) < 2Ce, and n i(*) < 2Ce(l + 7 In (1 + i)) 
Proof. Let T be the largest time such that 

(14.22) N 01 (t)= [ (1 + s)" 1 " 7 n i(s)fis + 1 < 2, for t<T 

Jo 

Then for t < T (|14.21|) holds and since, 

f'OO f'OO 

I (1 + s)^ 1 " 7 (l + 7 In (1 + s)) ds = 7 -1 / (1 + r) e" T oY = 27- 1 + 1 
Jo Jo 

it follows that 

N 01 {t) < 2Ce{2^- 1 + l) + 1 < 3/2, for t <T. 

Since A^i(t) is continuous this contradicts that T is the maximal number such that (|14.22l) holds. Thus 
T = 00 and (|14.21|) holds for all t < 00. □ 

This proves the first part of 1)14.6)1 and 1)14.7)) . The second part of 1)14.6)) follows from integrating 
the first using the boundary assumption ()14.2)) as in the proof of 1)14.3)1 . 

14.4 Proof of (jTOj) in case k = 1. 

We will now prove the first part of ()14.8)) for |/| = 1 assuming 1)14.1)1 - 1)14. 7)1 . 

Lemma 14.7. Suppose that the assumptions of Proposition hold and let F^y = F^ u 
be as in Lemma\3.°A Then 



(14.23) \ZF\ < Cet~ 1 {\dZh\ + \dh\) 

Proof. This follows from Lemma 113.31 □ 
Using the second part of Corollary I1U.3I 1)10.25)1 . and ()14.1)) - ()14.5)) and the previous lemma we get 
Lemma 14.8. If e > is sufficiently small then 

(14.24) (1 + t) ||3Z 7 /i(t,-)||L- < Ceil + tf £ (l + /* (1 + t)~ C£ ^ \\dZ T h(r, -)\\ L oo dr) . 

|j|<i - 70 |/|<i 

The estimate (pXH)) for |/| = 1 is now a consequence of the above lemma and the following technical 
result applied to m(t) = (1 + t) x - Ce Em<i \dZ I h{t, 01k 00 : 
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Lemma 14.9. Suppose that n\{t) > satisfies 

(14.25) m{t) < Ce(l + J (1 + r^niCr) dr^j 
Then 

(14.26) ni (t) <Ce(l + tf £ 
Proof. 

(14.27) Ni(t) = l+ [ (l + r) _1 ni(r)dr 

J o 

satisfies -/Vi(i) < Ce(l + r)~ 1 iVi(i). Multiplying by the integrating factor (1 + t) and integrating 
we get Nx(t) < iVi(0)(l + t) Ce = (1 + t) Ce and the lemma follows. □ 

14.5 Proof of in case jfe > 1. 

We will now use induction to prove the first part of (fl"4~8|) for 1 1\ = k + 1 assuming that (|14.1|) - (|14.5j) . 
the first part of ifTHfl) . ((T4~7|) and the first part of <fT4~g|) for |/| < k hold. 

Lemma 14.10. Suppose that the assumptions of Proposition [7^. i| hold and let F^ v = F^ u (h)(dh,dh) 
be as in Lemma VS .'A Then 

(14.28) \Z I F\<Cet- 1 ^ \dZ K h\ + C ^ \8Z J h\\8Z K h\ 

\K\<\I\ |J|+|K-|<|Z|,|J|<|K-|<|J| 

Proof. This follows from Lemma 113.31 using (|14.1|) - (jl4.7l) . □ 
By Corollary ITlil 

(14.29) \n g z I h\<\z I F\ + (i + ty l Y, E \z J H\\dz K h\ 

\K\<\I\, |J|+(|JC|-1)+<|I| 

+C(l+q)^ Y ( E \Z J H\ll+ E I^Ur+ E \Z J "H\)\dZ K h\ 

\K\<\I\ |J|+(|K|-l) + <j/j |J'|+(|iC[-l)+<|/|-l |J"|+(|K|-l)+<|7[-2 

where {\K\ - 1)+ = \K\ - 1, if \K\ > 1, and 0, if \K\ = 0. Using Lemma ITP1 we get 
(14.30) 

(1+q)- 1 Y \Z J H\ LL + \Z J 'H\ LT +\Z J "H\<^- t + £ sup \8Z J " H(t,y)\ 

\J\<k,\J'\<k-l,\J"\<k-2 \J"\<k-2 */ 2 <ls/l<2* 

We hence obtain 

(14.31) \U g Z I h\<Ce{l + t)~ l Y \dZ K h\+ ^ sup \dZ J H(t, y)\ \dZ K h\ 

\k\<\i\ iji+|ati<|/|-i t l 2 ^\y\^ 2t 

Then we have proven that 
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Lemma 14.11. Let 

(14.32) n k (t) = (1 + t) J2 WdZ'hit, •) \\ LOB . 

\i\<k 

Then for \I\ = k: 

(14.33) \U g Z J h\ < C(l + ty 2 (en k (t) + n k ^{tf) 
By the first part of Corollary I1U.3I (|1U.24|) . it therefore follows that: 

Lemma 14.12. 

(14.34) n k (t) <Ce + C [ (1 + r)" 1 (en k (r) + n fe _i(r) 2 ) 

J o 

Our inductive hypothesis is nfc_i(t) 2 <Ce 2 (l + t) c ' £ so the bound n k (t) <Ce(l + t) 2Ce follows from: 
Lemma 14.13. Suppose that 

(14.35) n k {t) < Ce(l + tf £ + Ce I (1 + r) _1 n fc (r) (ir 

•/ o 

then 

(14.36) n fe (t) < Ce{l + t) 2Ce . 

Proof. Let iV fc (t) = J *(l + r)- x n fc (r) dr. Then |JV fc (t)| < Ce(l + t)~ l {{\ + t) Ce + iV fe (t)). Multiplying 
by an integrating factor gives (N k (t)(l + ty 2Ce )' < Ce(l + 1)- 1 -^ so N k (t)(l + t)~ 2Ce < C and hence 
N k (t) < C(l + t) 2Ce and n fc (t) < 2Ce(l + t) 2Ce . □ 

This proves the first part of (|14.8[) . The second part of (|14.8f) follows from integrating the first and 
using the boundary assumption ()14.2)) as in the proof of 1)14. 3|) . 

15 Energy estimates for Einstein's equations 

Recall the definitions 

(15.1) E N (t) = sup V f IdZ 1 ^ 2 , 

o<r<t mN Jx t 

In this section we prove the following theorem. 

Theorem 15.1. Assume that g = h + m satisfies both Einstein's equations and the wave coordinate 
condition for < t < T. Suppose also that for some < 7 < 1/2 we have the following estimates for 
< t < T: 
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1. For all multi-indices I, \I\ < N/2 + 4 

(15.3) \dz ! h\ + (i + Mr 1 !^! + (i + t)(i + igr^az^i < ce(i + *)~ 1/2 ~ 7 (i + Mr 1/2_7 , 

S. For all multi-indices I, \I\ < N 

(15.4) |Z 7 tf(s,g,u;)| < Ce(l + for q = 1, 

3. 

(15.5) + (1 + M)" 1 ^™ + (1 + < Ce(l + t)" 1 , 

^. For a// multi-indices /, |/| < iV/2 + 2 

(15.6) |<9Z J /i| + (1 + < Ce(l + t)- l+Ce , 

5. 

(15.7) E N (0) < e 2 . 

Then there are positive constants C k independent ofT such that if e < C^ 2 we have the energy estimate 

(15.8) E k {t) + S k {t) < I6e 2 (l + tf" e , 
forO<t<T and for all k<N. 

Remark 15.2. Once again we recall that our estimates hold simultaneously for the tensors h and 
H = —h + 0(h 2 ). We shall freely interchange h and H in the proof below. 

Proof. Recall that the components of the tensor h^ u = g^ u — m^ v satisfy the following wave equations: 

(15.9) tf xP d a d p hp,=Fr,„ 

Fp, = P(d„h, d v h) + Q„ v {dh, dh) + G„ v {h)(dh, dh). 

where 

(15.10) P{d„K 8 u h) = ^m aa 'd^h Qa , mW'dvhpp - ^m aa ' m^' d^Kp d u h a ,p, 
We prove the desired estimate by induction on k. We first establish the estimate 

(15.11) E (t) + S (t) < 8e 2 (l + tf 0£ 

for some constant Co- After that we shall assume that the statement ()15.8j) for k < N' — 1 and prove 
the corresponding statement for k < N' with some constant CV'- We shall base our argument on 
the energy estimate (|11.3|) for the solution of the wave equation O g 4> = F proved in Proposition 1 1 1 . ll 
Observe that the conditions of our Proposition on the tensor h = g — m imply the assumptions of 
Proposition II 1 . ll for the metric g. 
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15.1 The case of N' = 0. 

In this section we prove the basic energy estimate for a solution of the equation ()15,9|) . 

□ a V = F i*> : = P{d^,d v h) + Q^(dh,dh) +G tw (h){dh,dh). 
Recall that according to (|13.10|) of Lemma 113.31 we have a pointwise bound 

\F\<\dh\ 2 TU + \8h\\dh\ + h\dh\ 2 
Using the assumptions of the proposition we infer that 

(15.13) |F|<e^i 
\ j 1 + i 

Therefore, the energy estimate (|15.12|) with = h^ v implies that 

Using the smallness assumption on the initial data and the Gronwall inequality this, in turn, leads to 
the desired estimate (fl5.11|) . 

E (t) + S (t) <8e 2 (l + t) Co£ 

15.2 The case of N' = 1. 

To facilitate the exposition we first consider the case N' = 1. We start by noting that according to 
(T7TH|) of Corollary [71H we have that 

^gZh^y = ZF^ U + D^, 
where the term = \H g Zh^ u — ZDgh^ v satisfies the estimate 

| D |<(^M + i^i r+ ^ ) E | 8 ^| 



|/|<l 



Recall that the tensor H a @ = —h al3 + 0(h 2 ). Thus using the assumptions on h of the proposition we 
derive that 

I ~ i + t 

\I\<1 

On the other hand, inequality ()13.11|) gives the estimate 

\ZF\ < (\dh\ru + \9h\ + \h\ \8h\) (\dZh\ + \dh\) + C\dh\ \dZh\ + C\dh\ 2 \Zh\ 
Using the assumptions of the proposition we conclude that 

{dZ 1 ^ \BZh\ 



\ZF\ = \(Z + Cz )F\<e^ l ^ T f+e 



Tiii L + t (l+t)3(l + |g|)a- 
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Now using the energy estimate (|15.12|) with <\> = Zh^ and F = ZF^ + we obtain 
t r „, r rt r \a7lu\2 rt 



1 f lOZ 1 ^ 2 [ l f \BZh\\dZh\ 



7o 7e t (i + ki) 1+27 - V M l + t Jo MtHi + \ q \) 



where we used the Cauchy-Schwarz inequality to pass to the last line. Combining this with the energy 
inequality 1)15. 14|) we infer that if Ce < 7/2 then 



|J|<1 * Ii1<l ''i'T \ X T I^M |I|<1 |/|<1 * ' 

The desired estimate 

E 1 (t) + S 1 {t) < 16e 2 (l+t) Cl£ 
now follows from the Gronwall inequality and the smallness assumption on the initial data. 

15.3 The case of TV' > 1. 

In what follows we assume that we have already shown that 

(15.16) E N ,_ x (t) + Siv-i(i) < 16e 2 (l + tf"'-^ £ , 
and prove that there exists a constant Cat/ such that 

(15.17) E N >(t)+S N >(t) < 16e 2 (l+t) c ^' £ , 
We start this section by writing the wave equation for the quantity Z 1 h liv with |/| = N' 



where 



\I\ g Z 1 h iliV — Z 1 F^y + D 1 ^, 



E)\iu — OgZ 1 h^ u — Z I \3gh l j iU 



We apply the energy estimate (|15.12j) with the functions (j) = Z 7 /i^ and F = Z 1 F^ + D 1 
(15.18) 

Note that we can estimate 
(15.19) 

( I {\Z I F\ + \D I \)\dZ I h\dxdt< ( ^—\dZ I h\ 2 dxdt+ [ [ e-^l+^d^Fl 2 + \D I \ 2 )dxdt 
Jo J Jo 1 + 1 Jo J 

Here the first term is of the type that appears already in the energy estimate (|15.18|) , Thus it remains 
to handle the second term. 
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According to (|7.21|) of Corollary 17.61 we have that 

|/| 

(15.20) D^^Dl, 

(15.21) D' t = Dl 1 + D' k 2 + D' k 3 + D' k \ 

(15-22) |^|<£ £ ttttm i^"!' 

\K\=k |J|+(|Jf|-l) + <|/| m 

(15-23) !^ 2 |<£ E T^l^l' 

|jf|=* |j|+(|jc|-i)+<[j| m 

(15-24) |^f|<E E TTkfl^l' 



|X|=fc |J|+(|Jf|-l)+<|/|-l 

1 + l?l 



(15-25) l^ 4 |<E E TZU\ dZKh \> 



\K\=k |J|+(]if|-l)+<|J|-2 

The estimates for Dl with k < N/2. We must now estimate 

(15.26) J j ' £- l (l + t)\D{\ 2 dxdt 

Since k = \K\ < N/2 in (jl5,22j) -( )15,25|) it follows from the assumptions in the theorem that we can 
estimate 

(15.27) e^il + t)\dZ K h\ 2 < mm f ' 



.(1 + i) 1 -^' (l + |g|) 1+2 ^ 
and it thus suffices to estimate 

(15.28) j'Js'\l + t)\D-^ ixd t < g f j + |g £ |)1+27 



\j\<\iy 



1 '■ \Z J H\ 2 



(15.29) JJ e-Hl + t)(\D^ + \Dir)dxdt< E J Q J (1 + t y-Ce + frjja dx dt > 



J|<|J|-i 

(15.30) [ / B -Hl + t)|i^| a cfactt< E /„ / min ( (1 + ly-cs » (1 + [g V +27 ) ^ 



„u . |j|<i-fr u " 

Lemma 15.3. Let f be a smooth function satisfying the condition 

(15.31) \f\<e(l + t)-\ for q = l 
Then 

(15.32) [ [ , g |N1 ^„ ^ , ,.- ctedt < / ^-rrrr- / l<9/f da; dt + e 

1 ; 7o J (i + k!) 1+27 (i + i + kl) 2 ~7o (i + 1+27 7 
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and 



i/i 2 r /■ ^i^/i 2 



mm 



(1 + tf-Cz ' (1 + [ g |)l+27 ; (1 + |g|)2 ~ y Q _/ (1 + | g |)l+2 7 



dt < / , „ , dxdt 



(15 ' 34) +E V„ (I+iF** 

Proof. We shall repeatedly use the Poincare inequality (|9,1|) of Lemma 19. II 

/• i/(^)i 2 ^ < /• i„2^ , r \d r m\ 2 dx 

(i5 - 35) y (t+mf^~ y m ds+ y 

which holds for any value of a > —1/2, a ^ 1/2. In particular, using (|15.31[) . we obtain that 

f \f\ 2 dx ^ o f \f\ 2 dx 
(15-36) / -Jl\-—< £ 2 + / 



(i + kl) 2+2CT ~ 7 (i + kl) 



2a 



The estimates (|15.32|) and (|15.33|) now follow from (|15.36|) with cr = 0. □ 
We now note the following generalization of (|15.35j) 

\f(x)\ 2 dx ^ E /' lrr » IC , [I'lfirf.Lr 

(i + \ q \) 2 ~ (i + ty-cz 

St -^(t+i) 



(15.37) /^( I _^, I _^)m^< I _| ra / + 



(i + M) 1+27 



The proof of (|15.37f) can be reduced to ()15,35|) by subtracting a term which picks up the boundary 
value. We define 

(15.38) / = /-/. where J(r,u) = f((t + l),u) X (r/t) 

and x( s ) = 1) when 3/4 < s < 3/2 and x( s ) = when s < 1/2 or s > 2. Then 
(15.39) 



mm 



\f(x)\ 2 dx < f \m£dx : f\f(x)\ 2 dx 



+ (i + [ g |)i+ar; (l + |g|) 2 ~ J (l + |g|) 3 + 2 -r {l + t) 1 -^ J (l + \q\) 2 

St St St 



We now apply (|15.35|) to the function /, which vanishes at r = t + 1, and observe that 



s< 
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On the other hand, 

|2 



which proves (|15.37j) . 

Using the lemma above with / = Z J H, together with ()15,28|) . (|15.29|) and the assumption that 
En'-i < 16(1 + t) c n'-i £ we see that we can estimate 



£ J e-\l + t){\D[}\ 2 + \D{ 3 \ 2 + \D[^) dx dt < jf {1 + £ t)1+2 , E N ,(t) dt + e 2 jf* 



—e dt 



ft 

dt 



-Ce 

JU \^ i "J 

for all < N/2. 

It thus remains the term ()15.3U|) containing D 12 . We shall use the version of the Poincare inequality 
(|15.34|) to create the term d q (Z J 'H)ll, which can be then converted to a tangential derivative of Z J H 
via the wave coordinate condition. However, in order to implement this strategy we modify the term 
Z J Hll according to Lemma ll2.2l We recall the notation 

(15.42) Hf lv = Z J H^+ Y, 4l zJH ^ 

\J'\<\J\ 

If \J\ < N' then the lower order terms in the right hand side of (|15.42j) may be estimated using (j!5.29f) 
and (|15.33j) as before. According to Lemma ll2.2l and the pointwise estimates in (|15.6|) and (j!5.4jl 

(15.43) \d r H J LL \ < Y \dZ J 'H\ + ^ \Z Jm H\ ■ ■ ■ \Z J2 H\ \dZ Jl H\ 

\J'\<\J\ \.h\+..+\,J m \<\,J\,m>2 

< Y \dZ J> H\+ Y \Z Jl H\\dZ j2 H\ 



\J'\<\J\ \Ji\+\M<W 

< 



V \07-rm i £(1 + |g|)1/2 " 7 |.^ffl i £lzJ ' Hl 



Hence 
(15.44) 

r J I 2 _/•*/•, e\dZ J 'H\ 2 e\dZ J ' H\ 2 e\Z J ' H\ 2 



+ TZ TTTTT- + TZ -rTTTr— ; TT77 dxdt 



+ M) 1+27 (i + t) 1+27 (i + t) 1+2 ^(i + |g|) 2 



If we use (jl5.33|) with Ce in the exponent replaced by 27 we see that the last term can be estimated 
by the second term from the right plus a term from the boundary: 



15 45) f £%HIl? dxdt - V /' l ( £{dZ ' J ' H? I ^ dZJ ' H? I ' c 2 ) dxdt 
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As we argued, when estimating ()15.3U|) we can replace \Z J H\n by the left hand side of Q15.42|) . After 
that we use the version of the Poincare inequality (|15.34[) applied to H[ L and this together with ()15.45j) 
gives 



(15.46) f + t )\ D k 2 \ 2 dxdt < eS N >(t) +eE N ,{t) + e 2 I 

Jo J Jo 

Summarizing, we have proven that 



(1 + 



l-Ce 



dt 



(15.47) J J e-\l + 1)\ D{\ 2 dxdt < eS N ,(t) + eE N ,{t) + s 2 J £ Ce dt, k < N/2 

This concludes the estimates in the case k < N/2. 

The commutator in case k > N/2. We isolate the case when \K\ = N' = \I\. We can estimate 
its contribution to the D J N , by the following expression: 

]D i ,< V ( \H\ + \ZH\ \ZH\ CC + \H\ CT \, K v \dZ*h\ 

where to pass to the last line we used pointwise estimates from (|15.5|) . (|15,3[) . and (|15.4j) . In the case 
when N/2 < k < \I\ we estimate the contribution of the corresponding term in Dj,, with the help of 
P5.6 1 as follows: 

E E E (^ft 

\K\<\I\ \J\<N/2 |y| |A'|<|/| V ' 

Therefore, 

(1 ,48, 1 1 e -(l + t W t? <.[J £ m J^L + | E | ^ 

Using the inductive assumption 1)15. 16|) we can therefore estimate 

(15.49) jTjT £ - X (l + t)|^| 2 dxdt < e jT ^1 rft + £ 2 J* + £ r ^_ 2C£ , iV/2<fc<iV' 
The inhomogeneous term. By (|1^.12|> 

(15.50) \Z : F\< {\dZ J h\ru\dZ K h\ru + \dZ J h\\dZ K h\) + ^ \dZ J h\ LT \dZ K h\ 

|j|+|if|<|/| \j\+\k\<\i\-i 

+ \dZ J h\\dZ K h\+ \Z Jm h\- ■ ■ \Z Ja h\ \dZ j2 h\\dZ Jl h\ 

\J\+\K\<\I\-2 |Ji|+...+|J m |<|/|,m>3 
The highest order terms with one of | J\, \K\ or equal to N = \I\ arc bounded by 

(15.51) {\dh\ ru + \Bh\ + \h\\dh\) ^ \dZ I h\ + \dh\ 2 ^ \Z I h\ + \dh\ ^ \dZ T h\ 

\I\=N \I\=N \!\= N 

.2 



" ^ l£v + U + ^(l + ^ ,£ + (1 + ^(1 + ^+7 lBZlHl 
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The remaining lower order terms are of the form 

(15.52) \dZ J h\\dZ K h\+ \dZ J h\\dZ L h\\Z K h\ 

\K\<N,\J\<N/2 \K\<N, \J\,\L\<N/2 

F F 2 
<7 r^rr- V \dZ K h\ + - — — --— V iZ'/ll. 



: Y \9Z K h\ + - —f- 

\K\<N v ' v |7|<AT 

It therefore follows that 



(15.53) 

/"* f 1/ wAJ„,2 , , / v- /"* f\dZ K h\ 2 \dZ K h\ 2 \Z K h\ 2 , , 

/. 1 8 (1 + t)|Z F| ~ ,g//o 7 -TT7- + (TTwW + (i + w4 + M? 



1/0 V ^ |/|<7V 



/i| (it 

|I|<~ 



. /"* /" |dZ A "/i| 2 \dZ K h\ 2 e 2 , , /•* ecft ^ , , l2l 

~ e I J + (TTW^ + (TT^ + J (i + r)^ ^ |az ^ 

Here, to estimate the last term in the first row we used ()15.33|) with — Ce in the exponent replaced by 
27, which produced a term similar to the first term of the first line plus a boundary term. Using the 
inductive assumption 1)15.16)) we thus obtain 

* E N i(t) dr 2 /** £dl ~ 



(1 + r) 



l-Ce 



(15.54) f f s-Xl + t^Ffdxdt^s [ ^l^L + eS N ,{t)+e 2 [ 

J0 JSr Jo 1 + 7" J 

The conclusion of the proof in case N' > 1 The inequalities (|15. 18(1 - 1)15. 19j) and 1)15 .47|) . 
1)15.49)1 and ()15.54)1 imply that for some constant C: 

(15.55) E N ,(t) + S N ,(t) < 8E N ,(0) + + + Ce f E ^ dT + C£ 2 /" gdr 

Jo 1 + 7" Jo (i + T ) 

If we now choose e so small that Ce < 1/9 we can move the second term on the right to the left and 
multiply by 9/8 to obtain for some new constants 

(15.56) E N ,{t) + S N ,(t) < 9JM0) + Ce f E »'^ dr + Ce 2 f — £ -% 

Jo l + T Jo i l + T ) 



Ce 



This can now be integrated using a Gronwall type of argument. If G(t) denotes the right hand side 
then we have 

Multiplying with the integrating factor we get 



: (G(t)(l 



+ tr C£ ) < 



dt\ y ' y ' J ~ 1 + t 



3 
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and hence if integrate and use that Celn(l + 1) < (1 + t) Ce , for t > (as is seen by differentiating 
both sides), and use that by assumption (|15.7|) G(0) < 9e 2 , we obtain 

G(t) < G(0)(1 + t) Ce + Ce 3 In (1 + t)(l + tf £ < 9e 2 (l + i) Ce + e 2 (l + i) 2Ce < 10e 2 (l + t) 2Ce 

Hence we have proven that 

E N ,(t) + S N ,(t) < We 2 (l + t) 2Ce 
This concludes the induction and the proof of the theorem. □ 

16 Geodesic completeness 

Having constructed a solution metric g = m + h of the Einstein equations we need to verify that the 
resulting space-time (M 4 ,<?) is causally geodesically complete. Let 

X(t) = (x°(r),x(r)) = (t(r),x(r) = (t(r),ru(r)) 

be a causal geodesic parameterized by the affine parameter r. Such geodesies satisfy the equations 

(16.1) X a (r) + T%[X{r))X^jri = 0, 

X(0) = Y, X(0)=£ 

where Y is the point of the origin of the geodesic X(t) and £ is the initial velocity satisfying the 
condition 

(16.2) g a p{Y)ef = ~A 2 <Q 
for some constant A. Condition Q16.2JI is preserved in time, i.e., 

(16.3) g aP (X(r))X a X^ = -A 2 

In the following lemma we show that a vector r] causal with respect to the metric g is "almost" causal 
with respect to the Minkowski metric m. 

Lemma 16.1. Let r/ be a causal 4-vector, i.e, 

(16.4) gaftrfif <-A 2 <Q 
for some non-negative constant A. Then 

(16.5) A+ \rf\ < 2|t7°|, Vi = l,..,3 
Proof. Expanding g = m + h we obtain from (|16.4|) that 

-I^ l 2 + EIV! 2 <N-(I^ I 2 + E^l 2 ) 

i=l i=l 
and the desired estimate follows provided that \h\ < 1/4. □ 
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We choose a future oriented initial velocity £, i.e., x°(0) > 0. 
Proposition 16.2. Assume that h = g — m satisfies the estimates 14 

\h\\dh\ + \dh\ TU + \dh\LL<et-\ 

\dh(t,x)\ < et^ 1 , for \x\<t/2 

Let X(t) is a future inextendible causal geodesic. Then the values of the affine parameter r span the 
interval [0, oo). 

Proof. We start by considering a time-like geodesic X(t). Reparameterizing, if necessary, we can 
assume that the constant A = 1 in (|16.3j) . Then equation (|16.3j) and inequality 1)16.5(1 with A = 1 
imply that for all r > 0. 

(16.6) £°(t) > - + \x{t)\ 

We removed the absolute value from x (r), since x°(0) > 0. This is the only part of the argument, 
which uses the fact that X(t) is a time- like geodesic. The case of a null geodesic will require an 
additional argument. 

Assume that X(r) is a time-like geodesic of finite length r*. We first observe that 

lim |-X"(t)| = oo 

t— »r* 

which means that X(t) escapes to infinity 15 in finite proper time t*. This easily follows from the stan- 
dard ODE theory. The inequality l|16.6|) implies that x (r) controls X(t). Thus to obtain contradiction 
it suffices to show that 

lim x°(t) < oo 

T— >T* 

Throughout this section we will use consistently use the notation x = t. We recall that 

r /3 7 = g a<T {dph la + d^hpa - d a hpj) 
Thus, expanding the metric g = m + h, 

x° - {2d p h 0l - doh^x^x 1 = h ■ dh ■ \X\ 2 

We further observe that 

(16.7) dph^x^x 1 = — (/io 7 x 7 ) — h^x" 1 

We now additionally recall that d q hLL is the only derivative of h that does not have the decay rate of 
at least (x ) -1 . Thus 

doh^xPX-y = dqhfrXPx-y + eO(t- 1 )\X\ 2 = d q hLL\X L \ 2 +eO((x°)- 1 )|A > | 2 

14 These assumptions are consistent with the decay estimates for h proved in Theorem 1 14. II 
15 viewed from the point of view of the global system of wave coordinates on R 4 . 
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The expression 



X± = X a L a = -x° + frx' = -4-(t - r) 

\X\ (XT 



Moreover, 

d q hLL = 4d q h 00 + eO((x )- 1 ) 
Furthermore, introduce £(x°/r) a cut-off function of the set r > x°/2. Then 

8 q h 00 = (1 - ()d q h 00 + Cd q h 00 = eO{r l ) + d q ((h 00 ) - (d q ((x°/r))h 00 

We compute 

d q ((x°/r) = (r" 1 +x°r- 2 )('(x°/r) < (x )' 1 

since r > x° /2 on the support of £'(x° /r). Thus d q hoo can be replaced by 9 g (C^-oo) at the expense of 
a term of order eO((x°) _1 ). Therefore, 



/-|2 i;i i, ,-.|2 . . m ,. ;, ,|,-.|2 , >, , .0^-1) |^|2 

L <7-L A o /at \-(rLiJ\rL i „/'-)/'" „,0\ — In 



<9 g /^|X^| 2 = H/toolgl 2 = H(C/ioo)!^| 2 + eO((x°)- 1 )| ( ? 



= 4— (C/ioo q) - <h 00 q - 4d L ((h 00 )X L X± - 4d n ((h 00 )X" X± + eO((a 
Here, 

/ l (JC(r))=%((r),t;(r), W (r)) 

where g = x° — r, v = x° + r, and u; = The advantage is that dnhoo, c^/ioo already decay faster 
than (x )- 1 and dn((x°/r) = 0, while \d L ((x°/r)\ < (x )- 1 . Thus 

d q hLL\X±\ 2 = ^-.{(hooq) ~ thooq + ^((x )- 1 ) 
It remains to analyze the term 

(16.8) q = -f-(x° ~ —x l ) =x° - x 1 — + r-Vlxl 2 - r" 2 |x • x| 2 ) 

dr r r 

From the geodesic equation 1)16.1(1 we can estimate 

\x a \ < \dh\\X\ 2 

Additionally, since on the support of C(x°/r), r > x°/2, we have that the last term in (|16,8jl multiplied 
by C^oo contributes at most 16 e(x°) _1 . Thus combining everything together we have 

|- (x° - 2/to 7 x 7 + (h 00 q) = 0{e{x°T l )\X\ 2 

We integrate this identity between proper times < r. Observe that \X\ < \x°\ and that 

(x°y 1 \x°\ = -^\nx 
dr 



3 This is the reason for introducing the cut-off function (. 



61 



Thus 



x°(r) <i°(0) + (2h ^ - Ch 00 q)\ T + e f In dr' 



It follows that 



x° 



|x°(0)|. 

Integrating one more time and assuming that x°(0) = i(0) = 1 we obtain that 

From this we conclude that the time x° = t remains finite with r. This concludes the proof for time-like 
geodesies. □ 

We now address the issue of null geodesies X(t), 

g a(3 X a X p = 

Examining the proof above leads to the conclusion that is suffices to establish that the condition 
x°(t) > is preserved in time. 

Lemma 16.3. For a future oriented inextendible null geodesic X(r) defined on the interval [0,7"*) we 
have x°(t) > for all r G [0, t*). 

Proof. Let tq < r* be the first time when x°(to) = 0. Fix a sufficiently small small constant c. Then 
there exists a small interval of size 5 such that 

0<±°(t)<c, VtG[t -5,to] 

and 

(16.9) ±°(t -6)=c 

We observe that (|16.5j) with A = implies that |-X"(r)| < 2|x°(t)| and therefore, 

\X(t)\<2c, VtG[to-5,t ] 
Integrating the geodesic equation (|16.1|) we obtain 

/"TO 

\x°(t ) -x°(t -S)\ < / \T\\X\ 2 <ec 2 5 

J TO— 6 

Thus, using ([roU|) . 

i°(r ) > c-ec 2 5 > 

Contradiction. □ 
This completes the proof of Proposition 116.21 

We have shown that all future inextendible causal geodesies X(t) exist for all values of the affine 
parameter r G [0, oo). This means that the constructed space-time is future causally geodesically 
complete. Next we establish that all future oriented causal geodesies escape to infinity. 
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Proposition 16.4. Let X(t) be a future oriented causal geodesic. Then 
(16.10) lim \X(t)\ = oo 

r^oo 

Proof. The inequality (|16.6|) immediately gives the desired result for time-like geodesies. Recall that 
by Lemma ll6.3l we have that x°(t) > and thus x (r) is monotonically increasing in r. We now argue 
by contradiction. Assume that for all r > 

|X(r)| <C 

for some potentially large constant C. Then there exists a time to such that 

to = lim x°(t) 



Set To be the value of the proper time r for which t(ro) = t^ — 8 for some small constant 5. Integrating 
the geodesic equation we obtain that for t > tq 

(16.11) x°(t) = x°(t )+ T \T\\x°\ 2 dT' <x (t )+e f x° dt' < x° (r ) + e5 sup i°(r) 

J TO Jto—8 T0<t'<T 

Thus for any r > tq 

(16.12) x°(t) < 2i°(r ) 

Choosing a sequence of times To — > co such that x°(to) — > (such a sequence must exist, otherwise 
x°(r) — > oo) we infer from ([16.12J1 that 

x°(t) -> 

asr-> oo. We can then choose small constant c, <5 such that i(ro) = to — ^ an d 

i°(r ) = C) x°( T ) < c 

for all r > tq. Returning to 1)16.11(1 we see that 

\x°(t) — c| < s5c 

Thus 

x°(t) > - 
\ ) - 2 

for all r > tq and we obtained contradiction. □ 
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